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We provide the analysis of charmless two-body B — > VP decays under the framework of the soft- 
collinear-effective-theory (SCET), where V(P) denotes a light vector (pseudoscalar) meson. Besides 
the leading power contributions, some power corrections (chiraly enhanced penguins) are also taken 

into account. Using the current available B —> PP and B —> VP experimental data on branching 
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fractions and CP asymmetry variables, we find two kinds of solutions in \ fit for the 16 non- 

perturbative inputs which are essential in the 87 B —> PP and B — > VP decay channels. Chiraly 

enhanced penguins can change several charming penguins sizably, since they share the same topology. 

However, most of the other non-perturbative inputs and predictions on branching ratios and CP 

asymmetries are not changed too much. With the two sets of inputs, we predict the branching 

fractions and CP asymmetries of other modes especially B a —> VP decays. The agreements and 

differences with results in QCD factorization and perturbative QCD approach are analyzed. We 

also study the time-dependent CP asymmetries in channels with CP eigenstates in the final states 

' and some other channels such as B°/B° — > 7r ± p =F and B°/B° — > K K*^ . In the perturbative QCD 

approach, the (S — P)(S + P) penguins in annihilation diagrams play an important role. Although 

they have the same topology with charming penguins in SCET, there are many differences between 
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<— j , the two objects in weak phases, magnitudes, strong phases and factorization properties. 

> ' I. INTRODUCTION 

m 
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Studies on B decays are mainly concentrated on the precise test of the standard model (SM) and the search for pos- 
sible new physics (NP) scenarios. To map out the apex in the unitarity triangle of the Cabibbo-Kobayashi-Maskawa 
■ (CKM) matrix, many precise experimental data together with reliable theoretical predictions are required. In charm- 
less two-body non-leptonic B decays, the main experimental obscrvables are branching ratios and CP asymmetries. 
To predict these observables, one has to compute the hadronic decay amplitudes (M\M2\Oi\B) , where Oi is typically 
a four-quark or a magnetic moment type operator. Since three hadronic states are involved in these decays, the 
predictions on these observables are often polluted by our poor knowledge of the non-perturbative QCD. Fortunately, 
it has been suggested that in the m& — > oo limit, decay amplitudes can be studied in a well-organized way: they can be 
factorized into the convolution of non-perturbative objects such as B to light form factors and decay constants of light 
pseudoscalars/ vectors with perturbative hard kernels. In recent years, great progresses have been made in studies 
of charmless two-body B decays. These decays were investigated in the so-called naive factorization approach 
and the generalized factorization approach 0, 0, [f], @, 0]. At present, there are three commonly-accepted theoretical 
approaches to investigate the dynamics of these decays, the QCD factorization (QCDF) 0, Q, the perturbative 
QCD (PQCD) H^EaEil, and the soft-collinear effective theory (SCET) Q,Q- Despite of many differences, all of 
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them are based on power expansions in Aqc_d/tti6, where mj, is the 6-quark mass and Aqcd is the typical hadronic 
scale. Factorization of the hadronic matrix elements is proved to hold in the leading power in Aqc d /frib in a number 
of decays. 

In the present work, we will focus on the SCET. The matching from QCD onto SCET is always performed in two 
stages. The fluctuations with ofT-shellness O(m^) is firstly integrated out and one results in the intermediate effective 
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theory. At final stage, we integrate out the hard-collinear modes with off-shcllncss O^^Kqcd) to derive SCET//. In 
B — > M1M2 decays, both of the final state mesons move very fast and are generated back-to-back in the rest frame 
of B meson. Correspondingly, there exist three typical scales: the b quark mass mb, the soft scale Aqcd set by the 



typical momentum of the light degrees of freedom in the heavy B meson, the intermediate scale y/rrib Aqcd which 
arise from the interaction between collinear particles and soft modes. SCET provides an elegant theoretical tool to 
separate the physics at different scales and factorization for B — > M\Mi was proved to hold to all orders in a s at 

leading power of I /nib Q H 0, EE 

|2fJ. After integrating out the fluctuations with off-shellness m|, one reaches 



the intermediate effective theory SCET/, in which the generic factorization formula for B — > M1M2 is written by: 

(MiM 3 |0<|B) = T{u) ® <j> Ml (u)C B ^ Ah + Tj{u, z) ® <j> Ml («) g> (f^ Ah (z), (1) 

where T and Tj are perturbatively calculable Wilson coefficients which depend on the Lorentz structure and flavor 
structure. Calculations for these hard kernel functions are approaching next-to-leading order accuracy 
22l . |23| |. In the second step, the fluctuations with typical off-shellness mbh-QCD are integrated out and one reaches 
SCET//. In SCET//, end-point singularities prohibit the factorization of Q, while the function Qj can be further 
factorized into the convolution of a hard kernel (jet function) with light-cone distribution amplitudes: 

C/0) =<f>M 2 (x)® J{z,x,k + )®<f> B {k+). (2) 

An essential question is whether power corrections in SCET can be analyzed in a similar way. It is almost an 
impossible task to include all power corrections, but we can include the relatively important one. Importance of 
chiraly enhanced penguins has been noted long time ago, and numerics show that chiraly enhanced penguins are 



comparable with the penguin contributions at leading power. Thus in both of QCDF P,S,G3] and PQCD 
approaches, it has been incorporated into the decay amplitudes besides the leading power penguins. In SCET, the 

1,Q. 
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complete operator basis and the corresponding factorization formulae for this term are recently derived in Ref. [2c 
A new factorization formula for chiraly enhanced penguin was proved to hold to all orders in a s , and more importantly 
the factorization formula does not suffer from the endpoint divergence. In the factorization formula, a new form factor 
named £ x and a twist-3 light-cone distribution amplitude (j) pp are introduced. 

In Ref. 25], one phcnomenological framework is introduced, in which the expansion at the intermediate scale 



fXhc = \/ m bh-QCD is not used. Instead the experimental data are used to fit the non-perturbative inputs. This 
method is very useful especially at tree level, since the function T{u) is a constant and Tj(u, z) is a function of only u. 
Thus only a few inputs are required in decay amplitudes. In this framework, an additional term from the intermediate 



charm quark loops, which is called charming penguin [20j, |25|, [26J, |27|, [28|, l29(, is also taken into account. Charming 



penguins are not factorized into the LCDAs and form factors, since the heavy charm quark pair can not be viewed 
as collinear quarks. They are also treated as non-perturbative inputs. This method is first applied to B — ► Ktt, 
B — > KK and B — > tttt decays 1251 . Subsequently, it is extended to charmless two-body B — > PP decays involving 
the iso-singlet mesons 77 and 77' [301 ] - 

In the present work, we extend this method to the B — ► VP decays. We will use the wealth of the experimental 
data to fit the non-perturbative inputs (in our analysis, we also take the B — > PP decays into account). In doing 
this, we would assume SU(3) symmetry for form factors and charming penguins to reduce the number of independent 
non-perturbative inputs: there are totally 16 non-perturbative inputs to be determined. Utilizing the meson matrices, 
we give the master equations for the hard kernels for B — > M\Mi decays. After analyzing the B — > VP decays at 
leading power, we take part of chiraly enhanced penguin into account. With the chiraly enhanced penguins taken 
into account, we find most of the 16 inputs are not changed sizably except charming penguins. Flavor singlet mesons 
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r\ and 77' receive additional contributions (gluonic contributions) from higher Fock state component. In Rcf. 3tJ, the 
gluonic form factors and gluonic charming penguins which are responsible for B — > PP decays are fitted using the 
related experimental data. Since there are not enough experimental results, the authors find two solutions for these 
inputs. This situation is changed when considering B — > VP decays since we have more data to give more stringent 
constraint. Incorporating the B — ► VP experimental results for branching fractions and CP asymmetries, we find 
that our results are consistent with their second solution. We find two solutions for the inputs only responsible for 



B — > VP decays. One of the solutions for B — > V form factors are smaller than those given in Ref. 23|, where the 
B — > plPl data (pl denotes a longitudinally polarized meson), B — > p°p~ and B — > p + p~ branching ratios and 
CP-asymmetries S p + p _ and C p + p -, are used. Our second solution for B — ► V form factors is more consistent with 
them. Generally speaking, charming penguins in SCET have the similar role with (5 — P)(S + P) annihilation penguin 
operators in PQCD approach. Both of them are essential to give the correct branching ratios in these two different 
approaches. But there are indeed some differences in predictions on other parameters such as direct CP asymmetries 
and mixing-induced CP asymmetries. We also make some comparisons between these two objects. 

The paper is organized as follows. B — > VP decay amplitudes at leading power are briefly given in Sec. HU What 
followed is the factorization analysis in which chiraly enhanced penguins are taken into account. In sectionUI] utilizing 
the rich experimental data on branching fractions and time-dependent CP asymmetry observables, we give two kinds 
of solutions for the 16 non-perturbative parameters responsible for B — > PP and B — ► VP decays at the leading power 
accuracy. With the inclusion of chiraly enhanced penguin, most parameters remain unchanged except the charming 
penguin parameters. Predictions on branching fractions and other observables, including direct CP asymmetries, 
time-dependent CP asymmetries and ratios of branching fractions, are given subsequently. A comparison between 
charming penguins in SCET and annihilation diagrams in PQCD approach is presented in Section[V] Sec. IVII contains 
our conclusions. In appendix[Aj we give the master equations for the hard kernels in both b — > d and b — > s transitions. 

II. B — > VP DECAY AMPLITUDES AT LEADING POWER IN SCET 

In this section, we briefly review the factorization analysis at the leading power and collect the corresponding leading 
order short-distance coefficients. The weak effective Hamiltonian which describes b — ► D [D = d, s) transitions are 



3l|: 



n *ff = ~k E v qb v q * D [c\ol + c 2 o q 2 -v tb v t * D [ J2 c >°> + H - C -' (3) 

V 1 ^ q=u,c i=3 ' 

where V^,(m are the CKM matrix elements and in the following we will also use products of the CKM matrix elements 
\ q J > (q=u,c,t) defined by A, = V qb V* f . Funct ions Oi (i = 1, 10, 77, 8g) are the local four-quark operators or the 
moment type operators: 

• current-current (tree) operators 

Of = (q a b a )v-A{Dpq p ) v _ A , O q 2 = {q a bp) v ^ A (D l3 q a ) v _ A , (4) 

• QCD penguin operators 

3 = (D a b a ) V - A '^2%q' )v-A, 4 = (Dpb a ) V -A^2(^ a ^)v-A, (5) 

q' q' 

5 = (D a b a ) V -A^2(q'0'f])v+A, O e = {Dpb a ) v _ A ^{q' a q' p ) v+A , (6) 
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• electro-weak penguin operators 

O-i = ^(-Dqfea)v-^y^eg/ (cf q0)v + A, 
<?' 

q' 

• magnetic moment operators 

Or^ = -^D a a^P R b a F^, 8g = -^D a a^P R T^b G%, (9) 

where a and f3 are color indices and q' are the active quarks at the scale mb, i.e. q' = (u,d, s,c,b). The m,b is the 
b quark mass and we use nib = 4.8 GcV. The left handed current is defined as (q' a q'p)v-A = q' a lv(l — J5)q'n an d 
the right handed current (q' a q'g)v+A — QaJvO- + l5)q'p- The projection operators arc defined as Pl = (1 — 7s)/2 
and Pr = (1 + 75)/2. The electro-weak penguin operators 09,10 can be eliminated using e q qq = uu + cc — ^qq- In 
the following, we will work to leading order in a s (mi>). In the naive dimensional regularization (NDR) scheme for 
otsijnz) = 0.119, a cm = 1/128, m t — 174.3 GeV, the Wilson coefficients Ci at leading logarithm order for tree and 
QCD penguin operators are 

Ci- 6 (m b ) = {1.110,-0.253,0.011,-0.026,0.008,-0.032}, (10) 

while the Wilson coefficients for electroweak penguin (EWP) operators are: 

CV-lo(m 6 ) = {0.09, 0.24, -10.3, 2.2} x 10" 3 , (11) 

and for the magnetic operators C7 7 (mf,) = —0.315, Cs g (mi,) — —0.149. We have used the sign convention for the 
electromagnetic and strong coupling constant as = — igT a A a ^ — ieQfA^, so that the Feynman rule for the 
vertex is igT a j^ + ieQfj^. 

In the present work, we will adopt the notations as in Ref. and use A = \J Aqcd/ith,. The emitted quark and 
anti-quark mainly move along the direction n+ and the recoiling meson is moving on the direction n_ , where n± arc 
two light-cone vectors: n± = and n + ■ n_ = 2. The matching from QCD onto SCET are always performed in two 
stages. We will first integrate out the fluctuations with off-shcllncss 0(ml) to give the intermediate effective theory. 
At final stage, we integrate out the hard-collinear modes with off-shellness 0(mi,AQCD) to derive SCET// 



3 _ ^ 

Os = ^(D(3b a ) V - A ^2e q '{q a q'p) 



V+A, 



°io = -^{Dpbodv-A^eq'iq'aq'^v-A, 
q' 



(7) 
(8) 



A. Matching onto SCETj 



To study the decay amplitudes of B — > MiM 2 decays in SCET, we first consider the possible operators using the 
building blocks. The power counting rule for these blocks has been given in Ref. [32j. Integrating out the hard scales 
with typical off-shellness m^, the electro- weak operators can match onto two kinds of operators in SCET where the 
situation is similar with that in B to light form factors: the first kind of operators involve four quark fields while 
the second one involves an additional transverse gluon field. For flavor-singlet mesons, one needs to consider the 
operators which are composed by two gluon fields. Then the leading power operators responsible for b — > s transitions 
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are chosen by: 



Qi°J(t) 



(sW c2 )(tn-)^(l- l5 )(W^u) 
{uW c2 )(tn_)*t-{l T ls){Wl 2 u) 
W2)(in-)^(1-7 5 )(W] 2 ?) 
(qWc2)(tn-)^-(lT K)(w! 2 q) 



[{uw cl )4+{i - l5 )h v ] 

[(sW cl )i+(l - l5 )h v ] 
[{qW cl )4+{\ - lb )h v ] 

[{sw cl )4 + {\~ lb )h v \ 



Q$(t) = m b ie^ u TT[[W^D^ c2 W c2 }(tn^)[w} 2 iDl c2 W c2 ]\ [{sW cl )4+(l - lb )h v ] 
with the trace over the color indices. The operators suppressed by A are given by: 



(12) 





1 

m b 


(sW c2 ){tn-)^—{l- lb )(wt 2 u) 




P±dW cl )(sn + )(l - 


■ is)K 


• 




1 

m& 


(uW c2 )(tn^^-(l T J 5 )(W} 2 u) 
ri—v 


'(sW cl )(W^i 


P±dW cl )(sn + )(l - 


■ ls)h v 


• 




1 

m b 


(sW c2 )(tn^)^—(l- l5 )(W} 2 q) 
ri-v 


\qW cl ){Wl x ii 


Dx c iW cl ){sn + ){l - 
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^(M) 


1 

m b 


(qW c2 )(tn_)^—(\Tls){w! 2 q) 


\sW cl ){Wt x il 


^LciW c i)(*n+)(l - 


lb)K 


1 




1 

m b 


\- S W c2 ){tn_)4^{l + l5 )(Wl 2 u) 


\uW cl )(W^ 


:Dl cl w cl )( sn+ )(i 


- 75) (sn + )h v 




1 

m b 


\sW c2 ){tn_)4^{l + l5 )(Wl 2 q) 


\qWoi)(W^t 


D^ cl W cl )( S n + )(l 


- l<o)K 


1 


Q^(M) 


= -2m h ie^ v Ti [[W^tD^ c2 W c2 ](tn^)[w} 2 


iDl c2 W c2 ] 







(sW cl )(W^ip^ cl W cl )(sn+)(l - l5 )h v 



(13) 



where the fields without position argument are at x = 0. The field products within the square brackets are color- 
singlet and we will neglect the colour-octet operators since they give vanishing matrix elements at leading order. The 
operators responsible for b — > d transitions could be directly obtained by replacing s quark fields by the corresponding 
d quark fields. Although the operators given in Eq. (fTB")) are suppressed by A compared with those in Eq. (fT2"|) . all of 
the operators in Eq. (fT2"|) and Eq. (fl"3"|) contribute to (MiAf2|0|£?) at the same power when matching onto SCET//. 
Hence the effective Hamiltonian are matched onto SCET/ by the following equation: 



H 



eff 



V2 I 



dt c,(<)O 4 (0) (t) + / dids bi(i, s)0^ ; (t, s) j> , 



(i), 



(14) 



with s = n + ■ p's = nigs, i = ra_ ■ qt = mst (p' and 5 are the momentum of the recoiling and emitted meson, 
respectively). We usually evaluate the Wilson coefficients c 2 (u) and bi(u,z) in momentum space which is related to 
the ones in coordinated space by: 



Ci(u) = I die 



(i), h(u,z)= die- im ^ ut+z %(i,s)- 



(15) 
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The tree level matching coefficients for the four-body operators in eq. (|12[) arc given by: 



„(/) 

'1,2 
„(/) 



„(/) 
-4,5 

„(/) 



4/) 



Ci } 2 + —^2,1 



-^rCg^io + ClQ,9 



,(/) 



-a; 



(/> 



-a; 



(/) 



r 1 

c 5 



C3,4 + C4 3 



1 



N~ c Ce ~l C7 ~2N c C 8 



The tree level matching of five-body operators leads to: 



°4,f 

°7 
°8 

6^ 



! L V uj 2 ' N c 
C 4 , 3 + ^(1-^)C73 

iVc V W3 / 

c 5 



.(/) 



3 r 



c 



10,9 



^3 ' 



1 

A^ 



C03 

U>2 > 



+ A 



2 L A'c V W3 / 



+ A^ 



(/) 



C 7 



W2 



3 1 rm b _ m b \ 
2 7 Ar c VcJ 2 



c 5 - ic 7 



AT C V w 2 W3 
o; s (m b ) / f _ 1\ r 2 + z 
' 16Cf u)V\ — z 



2 1 



A 7 ?/ (1-zu)(1-zm) 



(16) 



(17) 



where W2 = urns and W3 = —urns with w is the momentum fraction of the positive quark in the emitted meson, wig is 



the £?-meson mass. Cf = (A 7 ,? — l)/2Af c and AT C = 3. The one-loop corrections are given in Rcfs. 
The coefficients c| and b£ are zero at ©(c^), thus they are not relevant for the present study in which we concentrate 
on the leading order analysis. 

In SCET/, the matrix elements of of*' 1 ^ can be decomposed into some simple and universal ones defined as follows: 

due iut (/) Ml (u), 



(18) 



(M 1 |(x^c 2 )(in-)^(1 - 75)(^ C U)|0) = 
(Af 2 |T[(x^ c i)^+(l - 7S)^]|^) = m B (, 

(M 2 \T{(uW cl )(W f cl ifl ±cl W cl )(sn + )(l - l5 )h v ]\B) = -m% J dze imB *- a {j(z), 
where M<x is an arbitrary pseudo-scalar meson or vector meson except r\ and rf . 



B. Matching to SCETn 

The matching of SCETi onto SCETn is performed by integrating out the degrees of freedom with p 2 ~ Am&. To do 
so, it is useful to perform a redefintion of collinear fields: q — > Y s q, where Y s is a soft Wilson line. The SCET Lagrangian 
contains no leading order interactions between the collinear-2 and collinear- 1 fields after decoupling soft-gluons from 
collinear-2 sector by a field re-definition. Although soft Wilson lines still appear in the effective electro- weak operators, 
the Wilson line only appear in the combination of Y s h v . Thus the two kinds of collinear sectors decouple and the 
decay amplitudes factorize. 
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In SCET//, the end-point singularity prevents the factorization of C while the form factor Cj M ( z ) can be further 
factorized into convolution of light-cone-distribution amplitudes (LCDAs) and jet functions: 

CfW = fdk + dx(b+(k + )J(z,x,k + )^ M (x). (19) 

m B J 

At the lowest order, J(z,x,k + ) = 5(z — x)a s irCF /(N c xk+). 



C. Decay amplitudes involving flavor-singlet mesons r\ and r/ 
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34, M. In this 



For iso-singlet mesons r\ and 77', we adopt the Feldmann-Kroll-Stcch (FKS) mixing scheme 
scheme, an arbitrary iso-singlet biquark operator O can be written as a linear combination of O q ~ (uu + dd)/\/2 and 
O s ~ ss operators with the well defined flavor structure. Matrix elements of O — c q O q + c s O s between 77, 7/ states 
and the vacuum state can be parameterized by 

(0|O|tj) = c g cos <t> q (O g }-c s sin <j> s {O s ), (20) 
(0|O|V) = c q sin <j> q (O q )+c s cos <f>,{O a ), (21) 

where the four matrix elements (0|O 9iS |?7(')) are expressed by the two angles 4>q,s and two reduced matrix elements 
(O q>s ). Phenomenologically, one can neglect the OZI suppressed matrix elements and obtain cj) q — <fi s = 8. Thus, the 
mass eigenstates 77, 77' are related to the flavor basis through: 

77 = rj q cos 8 — r] s sin 8, 

rf = rj q sin 8 + r] s cos 8. (22) 

For these iso-singlet mesons r\ q and r\ Sl we need in addition more theoretical inputs which arise from the higher Fock 
state component: 

3' 
T 



ze^(7, ? (p)|Tr[[W c V^ c2 ^ c2 ](tn_)[Fy c V^ c2 ^ c2 ]]|0) = / due«*-y/C F ^ -f Vt & P (u), 



7^<77 s ( JJ )|Tr[[VFj 2i ^ c2 ^ c2 ](tn_)[^J 2 7Cl c2 ^ c2 ]]|0) = J due mt C F ^f^ P {u\ 

(( Vq \T[(xW cl )4 + (l - l 5 )h v ]\B)) g = V2m B ( g , 
(( Vs \T[(xW cl )4+(l - J 5 )h v ]\B)) g = m B ( g , 

{{i lq \T[{ X W cl ){Wl l iIp^ 1 W cl ){ S n + ){l -j 5 )h v ]\B)) g = -V2m| J dze™*™ '{j g (z), 
((i h \T[(xW cl )(W f cl tp ±cl W cl )( S n+)(l - lh )h v ]\B)) g = -m% { dze ml * z - s Cj g (z), (23) 



where only the gluonic contributions to B — > rj qi ij s form factors are shown. Please note that, our convention is 



BM 2 



30l | . where the form factors £ g and Qj g arc incorporated in the definition of 



different from the one used in Ref. 
Here we have separated them out and the two functions Cfy) 2 do not contain contributions from the gluonic term. 
This convention is more convenient when extracting the hard kernels using master equations given in the appendix. 

In SCET77, C g can not be factorized either for the presence of end-point singularity but Cj g i z ) is given in terms 
of the jet functions by: 



i-bm 1 j\ _ Mm 1 , fCj f .....+ 



G7(z) = ^-aH- dk + dxcb+(k + )J g (z,x,k + )^ M (x), (24) 



m B 4 v 3 J 

At the lowest order, J g (z, x, k + ) = 5(z — x)a s 2-K / '(N c k + ). 



8 



D. A summary of the factorization formulae 

In summary, the b — > s(d) decay amplitudes at leading power in SCET can be expressed by: 

A(B — > M±M 2 ) = -^m|{/ Ml Jducp M MTi(u)t BM * + f Ml Jdu0 Ml (u) JdzTuiu, z)C^ M2 {z) 
+f Ml J du4> Ml {u)T lg {u)C, BM * + f Ml J du<j>M, («) J dzT 1Jg (u, z)(f g Ah (z) 
+fl h J du& Ml (u)T?(u)( BM * + f Ml J du& Ml (u) J dzTfjiu, z)Cf A h (z) 
+fl h J du& Ml (u)T? g (u)( BM > + fl h J du& Ml (u) J dzT? Jg (u, z)Q B g M * (z) 
+\V'>A% lM * + (1 «- 2)}, (25) 

where A^, lMz denotes the non-perturbative charming penguins. Tj arc hard kernels which can be calculated using 
perturbation theory. In the appendix ^ based on the flavor structure of the four-body operators and five-body 
operators, we give the master equations for hard kernels Tj which utilize the coefficients given in Eq. (|16|) and 
Eq. p7|) . For distinct decay channels, one can easily evaluate the equation to obtain the corresponding hard kernels. 

In SCET, the factorization formula for B — > M\Mi is easily proved to hold to all order in a s : the amplitudes given 
in Eq. (|2"S")) have the form of a convolution of the universal light-cone distribution amplitudes and the perturbative 



hard kernels. Utilizing the perturbative expansion in a s (v"^A) for the jet functions and in a s {rrib) for the Wilson 
coefficients, one can predict the branching ratios, CP asymmetries and other observables for B — > M\Mi decays. One 
can also use another parallel method: the non-perturbative parameters can be fitted by experimental measurements 
on the B — > M\M% decays. This approach is especially useful at leading order in a s , since then the hard kernels T\{u) 
are constants, while T\j(u,z) are functions of u only. Furthermore, at this order terms with hard kernels Tfj{u,z), 
Tf(u), Tfj (u, z), Tf g (u) do not contribute at all. Thus the decay amplitudes of B — * M\M<2, decays at LO in a s ( m i>) 
are written by: 



A(B MiM 2 ) = ^|m| (f Ml Cj M2 J du4> Ml («)Tij(«) + (% Ah J duct> Ml (u)T 1Jg (u) 

+f Ml (7K SM2 + T la C BAh - ) + \PA% m * + (1 <-> 2) } , (26) 
where the four functions ( BMl , Q and 

= J dzCf M Hz), OT 2 = / dzC?*(*), (27) 

arc treated as non-perturbative parameters to be fitted from experiment measurements. 

In order to reduce the independent inputs, one can utilize the SU(3) symmetry for B to light form factors and 
charming penguins. In the exact SU(3) limit, only two form factors are needed for B — > PP decays without iso-singlet 
mesons: 

rBP _ /-Bit _ j-BK _ /-B S K /r, Q \ 

C(j) = C(j) - C(j) - C ( j) ■ (28) 

Besides these two form factors, there are two additional new non-perturbative functions Oj)g m decays involving 
iso-singlet mesons r\ q and r\ s . They are contributions from the intrinsic gluons. The B — -> V form factors are rather 
simple, since there is no gluonic contribution at all. The flavor SU(3) symmetry implies the relation for B — > V form 
factors: 

/-BV /-Bp j-BK* /-Bu) ^B,K" /-B s 4> (on\ 
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If the SU(3) symmetry is assumed for charming penguins, there are totally five complex charming penguins which 
depends on the spin and isospin properties of the emitted mesons and recoiling mesons: A^f , A^f , A^f , A^,A^. 
Acc lM2 denotes the charming penguins in which the M\ meson is emitted and the M 2 meson is recoiled. The two 
charming penguins A^ g , A^ only contributes to decays in which a iso-singlet meson is recoiled. 

With the assumption of flavor SU(3) symmetry for B to light form factors and charming penguin terms, the non- 
perturbative, totally 16 real inputs responsible for B — > PP and B — > VP decays are summarized in the following: 



c 



BP 



SJ ! 



Cg, C.Jg, C 



BV 



/■BV 
W > 



.4 



pp 



A 



pv 



A 



VP 



A 



pp 

ccg ) 



.4 



VP 
ccg ■ 



(30) 



III. CHIRALY ENHANCED PENGUINS 



Power corrections are expected to be suppressed by at least the factor Aqc^/to;,, but chiraly enhanced penguins 
are large enough to compete with the leading power QCD penguins as the suppression factor becomes 2/xp/mb, where 
lip ~ 2 GcV is the chiral scale parameter. Thus in both of QCDF and PQCD 0, Q Q approaches, it has 

been incorporated in the phenomenological analysis. In the framework of SCET, the complete operator basis and the 
corresponding factorization formulae for the chiraly enhanced penguin are recently derived in Ref. 23j, l24j and the 



amplitudes do not suffer from additional endpoint singularities. The factorization formula will introduce a new form 
factor C x and a new light-cone distribution amplitude 



As discussed in Ref. [23j, there are three different kinds of chiraly enhanced penguin operators in SCET/: Q 1 ^, 



Q£ x) and Q {2x) . 



The basis for the Q^'-type operators is given by: 



{lx) 

^ 1(9/9) 



i 



m b 



[(qW c i)(l - j 5 )h v 



(sW c2 )(tn-)^—i$ x (l + 7 5 )(W«k) 



,(ix) £ p 
'1(9/9)2 q ' 



(31) 



2(9/9) 

12 w iH contribute to B — > PP, VP, VlVl decays (here Vl denotes a longitudinally polarized 
There are in addition several operators omitted here, as they can only contribute to B — > VtVt 



These two operators Qx 2 ' w ^ contribute to B 



vector meson) 

decays (Vr denotes a transversely polarized vector meson). The second kinds of operators which are responsible for 
B — ► PP, VP, VlVl decays are given by: 

1 



{2x) 
^1(9/9) 



o {2x) 



Q 



(2X) 
3(9/9) 



-1 

m b 
x 

m b 
x 

ZT2 



(qWcl)- 



n + ■ id 

(sW c2 )(tn.)^(l- l5 )(W^q) 
1 



id ± ■ (W f cl iD ±cl W cl ){sn+){l+j 5 )h v 



(sW c 



cl 



id 



id±_ ■ (W^iD Xcl W cl )(sn+)(l + l5 )h v 



(uWc2)(tn^(l+ l5 )(W^ 2 u) 
(qW cl )(W 1 cl ip ±cl W cl )(sn+)(l - 75 )h v 

(sW c2 )(tn-)^ijF) ± (l + -y 5 )(W^q) 



(32) 



(33) 



Q 



(2X) 
4(9/9) 



-e {2x) 



plus operators with the same Dirac structure but different flavors, Qx$f u \ and Q^f uu y 
included, we have two additional operators Qi(j qq ) and Q^/gq) - ' 



(34) 
(35) 

If n_-iso-singlet operators are 
Operators contribute to B -> PP,VP,V L V L 
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decays, while operators which only contribute to B — > VpVp decays are also given in Ref. 23[ but omitted here, since 
we mainly concentrate on B — > PP and B — > VP decays. 

Matching from QCD to SCET/, one obtains the effective Hamiltonian expressed by the (lx) and (2x)-type operators 
contributing to B — > PP, VP, VlVl decays: 



H 



eff 



V2 



dt< 



$)Q [ i(F) (*) + / dtdsbf (F) (t, s)Q^> (t, a) 



2X), 



(36) 



where the indices run over the operator number i and possibilities for the flavors F for the Q^p)- <%(f) an< ^ "i(F) are 
the short-distance Wilson coefficients in coordinate space. At tree level, the corresponding coefficients in momentum 
space are: 



h x 



b l(ufu) 



b x 



b x 



Cio 



C 7 , 1 



a; 



(/) 



1 + UZ C 3 _ Cg_ 

uz [ N r 2N r ' % 2 



. ^ 2(fuu) =3X ( t n [C 7 



2(1 + uz) ( C 2 



3C< 



-^A1/) + ^aF ) )-(2C 1 A(/)-3C 10 aF ) ) 



2N r 



(fh 



N r UzN r 



(37) 



Matrix elements for these operators can be parametrized into the following universal distributions: 



(M(p)| 



(qW cl )^—id±_ ■ (W^iD ±cl W cl ){sn+)(l + l5 )h v 
n ■ id 



(sW c2 )(tn-)^—ip x (l + j 5 )(W} 2 q) 



|o) = - 



B) 

ifMUM 



6 



dze imBZ - s C^ M (z), 



due lut q>Z(u) 



(38) 



where [iu is the chiral scale parameter which is set to zero for vector mesons. Using equation of motion, the pseudo- 
scalar's light-cone distribution amplitude (ffp(u) can be related to ones defined in QCD 2J, |36j: 



4> P p{u) = 'iu 



2/3P f dv 



-<j>3p(u - V,U) 



(39) 



6 fpfJ-p 

In the Wandzura-Wilczek approximation, <p 3 p vanishes and one gets 4> P p(u) = 6m(1 — u) for the asymptotic form 
With the above matrix elements, generic decay amplitudes from the chiral enhanced penguin could be written as: 

MMi /mi 



A X (B — > MiM 2 ) = 

imp 

_/WMi_ 

3ms 

A* Mi fati 

3m p 
6m r 



3mp 



ducj>^{u)T x (u)C BM * 

BM 2 



dudzcj)^{u)T*{u,z)tf M *{z) 
du^(u)T x (u)^- 



dudz^(u)T x Jg (u,z)Cf g M2 (z) 



dudz4> Ml (u)T x {u,z)C BM2 {z) + (1 



4 



(40) 



where Cx( z ) can ^e expressed as convolutions of LCDAs and jet functions: 

CW = Jo *c JT dk+ i^5l^+(fc+)^( x ). 



(41) 



11 



Here J±(z, x, k + ) — S(x — z)ita s CF / '(N c xk+) at lowest order. 

As emphasized in section [Til the leading power SCET phenomenological analysis is very useful especially at tree 
level. It does simplify the analysis. Even taking into account the first four terms in Eq. (|40j) . the scheme for 
phenomenological studies will remain. But considering the chiraly enhanced penguins, the factorization formulae 
involves a new form factor £ x which can not be simplified into a normalization constant even at tree level. As shown 
in Rcf. 23[, the fifth term proportional to £ x is small which docs not give sizable contributions. Thus in our analysis, 



we neglect it and only consider the first four terms: 



BM 2 
Ug^Jg 



+ (1^2)}. (42) 

For B — > PP decays, the chiraly enhanced penguin takes a plus sign; while in B — > VP decays, when emitting a 
pseudoscalar meson, the amplitude take a minus sign; when a vector meson emitted, there is no contribution from 
chiraly enhanced penguin since [iy = 0. 



IV. NUMERICAL ANALYSIS OF B — > VP DECAYS 



A. Input parameters 



In the factorization formulae, we will use the following values for decay constants of the light pscudo-scalars and 
vector mesons ( in units of GeV): 



U = 0-131, f K = 0.160, f, lq = 1.07/* = 0.140, /„. = 1.34/* = 0.176, 
f p = 0.209, f K . = 0.217, U = 0.195, U = 0.231. 



(43) 



33 



34 



3J: 



The mixing angle between r\ q and rj s is chosen as 9 = 39.3 
angles, we use the updated global fit results from CKMfitter group 

V ud = 0.97400, V us = 0.22653, \V ub \ = (3.57±g:J?) x 10" 3 
V cd = -0.22638, V cs = 0.97316, V cb = (40.5i^) x 10~ 3 , 
\Vtd\ = (8.68±g;i|) x HT 3 , \V ts \ = (40.7±g-J) x 10" 3 , V tb = 
= pLTjffiT, 7 = (67.6^)°, e = (1.054±8lgg)°. 



35J . For the CKM matrix elements and CKM 



(44) 



For the inverse moments of light-cone distribution amplitudes for pseudo-scalar mesons, we use the same value as in 
Ref. 



30J: 



1q 



= 3.3, 



(x- 1 ) K = 3.2A, (x~ l ] 



K 



3.42, 



(45) 



where the inverse moment of vector mesons' light-cone distribution amplitudes are obtained utilizing the Gcgcnbaucr 
moments evaluated in QCD sum rules 



(x- l ) p = (!-% = 3.45, (x- 1 )j > = 3.54, (x' 1 ) K , = 2.79, (x" 1 ' 



K 



= 3.81. 



(46) 



For the chiral scale parameters, we use a universal value \xp = 2.0 GeV for pseudo-scalars and /Lty = for vectors. 
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The experimental data of B — > PP and B — > VP branching ratios, the direct CP asymmetries and the parameters 



in B°/B° 



decays (which are defined in Eq. (|73"1). Eq. (|74| and Eq. |75]1) are given by Heavy-Flavor- Averaging- 



Group (HFAG) [39( and Particle-Data-Group (PDG) [4Cf ■ The following mixing- induced CP asymmetries in B 
and B — > VP decays are also used in our analysis: 

-rifS(Ksr)') = 0.61 ± 0.07, -7] f S(K S TT ) = 0.38 ± 0.19, S^tt+O = -0.61 ± 0.08, 
-ri f S{<pK s ) = 0.39 ± 0.17, S{ir°p°) = 0.12 ± 0.38, 



PP 



where rjf is th 
since the expe 



-rifS(p K s ) = O.eitoil ± 0.09 ± 0.08 = 0.61 



d(s) 



+0.25 
-0.27' 



d{s) 



-r) f S(wK s ) = 0.48 ± 0.24, (47) 
is not used in this fitting, 



-a r> 



4 



d(s) 



d(s) 



b 



■U Ur 



<7 




FIG. 1: Feynman diagrams for chiraly enhanced penguins (left) and charming penguins (right). The two diagrams in the lower 
line only contribute to decays involving 77 or rj , where q — q . 



With these data for branching fractions and CP asymmetries, \ 2 fit method is used to determine the non- 
perturbative inputs: form factors and charming penguins. Straightforwardly, we obtain the two solutions for nu- 
merical results of the 16 non-perturbative inputs. At leading order and leading power accuracy, the first solution is 
(the charming penguins are given in units of GeV): 





= (12.8 ± 1.2) x 10" 


-2 


Cj = (7.2 ±0.7) x 10~ 2 , 


c v 


= (12.4 ± 1.8) x 10" 


-2 


Cj = (10.8 ± 1.9) x 10~ 2 , 


c 9 


= (-5.3 ±2.2) x 10 


-2 


Cj s = (-2.3 ±2.9) x 10~ 2 


\A P J\ 


= (48.1 ±0.6) x 10" 


-4 

3 


arg[A p c p ] = (167.5 ± 2.5)°, 


\A V J\ 


= (40.6 ±0.9) x 1(T 


-4 

; 


arg[Al P ] = (10.7 ±4.3)°, 


\A P J\ 


= (30.7 ± 1.3) x 10" 


-4 

; 


arg[A™] = (194.3 ± 4.6)°, 


\Accg\ 


= (38.4 ± 1.9) x 10" 


-4 

1 


arg[A p f g ] = (83.0 ±3.8)°, 


\A v J g \ 


= (23.0 ±2.4) x 10" 


-4 


arg[A v J g ] = (38.4 ±23.0)°, 
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and one can obtain the predictions for B — > P (here P denotes a pseudoscalar except r\ and rf) and B — > V form 
factors at tree level: 

F B ^ P = C P + Cj = 0-201 ±0.015, A B ^ V = C y + Cj = 0.232 ±0.037. (49) 

In the above equations (and also in the following), the uncertainties are obtained through the x 2_n t program. After 
including the chiraly enhanced penguin, the numerical results for these inputs are (the charming penguins are given 
in units of GeV): 





= (13.7 ±0.8) x 10" 


-2 


Cj = (6.9 ±0.7) x 10~ 2 , 


c v 


= (11.7 ± 1.0) x 10" 


-2 

5 


Cj = (11.6 ±0.9) x 10~ 2 , 


c 9 


= (-4.9 ±2.4) x 10 


-2 


C/ 3 = (-2.7 ±3.2) x 10~ 2 , 




= (40.0 ±0.6) x 10" 


-4 


arg[A pp ] = (165.2 ± 2.8)°, 




= (41.0 ±0.9) x 10" 


-4 


arg[Al P } = (11.9 ±4.2)°, 


A P J\ 


= (39.9 ± 1.0) x 10" 


-4 


arg[A p c v ] = (191.5 ± 3.6)°, 


A pp \ 

ccg \ 


= (37.7 ± 1.8) x 10" 


-4 

) 


arg[A pp ] = (88.3 ±4.1)°, 


A vp \ 

ccg \ 


= (25.3 ± 2.3) x 10" 


-4 

i 


arg[A v c p ] = (-18.7 ±12.3)° 



(50) 

which gives the predictions for B — > P and B — > V form factors at tree level: 

F B ^ P = 0.206 ± 0.004, A^ v = 0.233 ±0.017. (51) 

As shown in Fig. [T] chiraly enhanced penguins have the same topology with the charming penguins. The former 
two diagrams do not only contribute to decays without iso-singlet mesons rj or rf but also decays with these mesons. 
The two diagrams in the lower line only contribute to decays involving r\ or 77', where q = q' . The inclusion of chirally 
enhanced penguin will mainly change the size of three charming penguins A pp , A pp g , A p ff . Predictions for branching 
fractions and CP asymmetries will not be changed sizably. After including the chiraly enhanced penguins, the total 
X 2 /d.o.f for obscrvablcs B — > PP and B — > VP is 301/(86 — 16). If only the 55 observables in B — > VP decays are 
concerned, the total x 2 is 112. 

Besides the above results, there is another solution at leading power: 





= (13.4 ±0.3) x 10" 


-2 




= (5.8 ±0.4) x 10~ 2 , 


c y 


= (22.9 ± 1.3) x 10" 


-2 




= (6.6 ± 1.4) x 10~ 2 , 


c 9 


= (-10.3 ± 1.2) x 10~ 2 . 


( Jg - (5.8 ±1.5) x nr 2 , 


A-cc\ 


= (48.4 ±0.4) x 10" 


-4 


arg 


■A pp ] = (167.1 ±2.6)°, 




= (29.7 ±0.8) x 10" 


-4 

5 


arg 


\A^ C P ] = (159.3 ±6.9)°, 




= (44.9 ± 1.1) x 10" 


-4 


arg 


\Kc\ = (-10.5 ±2.9)°, 


A ccg\ 


= (38.4 ±2.2) x 10" 


-4 

5 


arg 


\A PP ] = (83.8 ±4.5)°, 


A vp \ 

ccg \ 


= (18.6 ±2.3) x 10" 


-4 

5 


arg 


\A V C P \ = (220.6 ±10.7)° 



(52) 



which gives: 



?B ^ P - 0.192 ±0.005, A B ^ V = 0.295 ±0.009. (53) 
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With the inclusion of chirally enhanced penguin, these inputs become: 



C P = (14.1 ±0.8) x 1CT 2 , 
C y = (22.7 ±1.7) x 1CT 2 , 
( g = (-10.0 ±0.9) x 10~ 2 . 
\A PP \ = (40.6 ±0.6) x 10" 
\A^ C P \ = (29.4 ±0.8) x 10" 
\A P ^\ = (33.5 ± 1.1) x 10" 
(37.8 ± 1.3) x 10" 



\A PP 



\A 



VP\ 



(18.3 ±2.4) x 10 



(j = (5.6 ±0.7) x 10" 2 , 
Cj = (6.5 ± 1.8) x 10~ 2 , 
Ci g = (5.1 ±1.1) x 10- 2 , 
arg[A pp ] = (164.9 ±2.8)' 
arg[A^ L p ] = (158.4 ±5.8) 
arg[A P y} = (-U.3±3.8] 
arg{A pp } = (87.5 ±2.1)°. 



arg[A 



VP ] 

ccg \ 



(225.6 ± 10.0) c 



(54) 



with the form factors: 



F 



0.198 ±0.003, A 



0.291 ±0.011. 



(55) 



The corresponding \ 2 = 271/(86 — 16) (x 2 for the 55 observables in all B — > VP decays is 69). Comparing the results 
in the leading order analysis and those with chiraly enhanced penguins, we can see that the charming penguins A pp 
and A P Y are changed sizably. It is reasonable since chiraly enhanced penguins and charming penguins have the same 
topology. The phase of AY C P is also changed sizably. It implies that the total statistical significance \ 2 1S n °t very 
sensitive to arg[AY c P ]. The large error in this parameter also confirms this feature. 

Using the two solutions for these non-perturbative inputs, we obtain two different kinds of predictions (labeled as 
This work 1 and This work 2) on branching fractions and CP asymmetries, where the chiraly enhanced penguins 
are taken into account. As we have shown in the above, the leading power results are not very different from these 
results, as the inclusion of chiraly enhanced penguins only amounts to a redefinition of charming penguins. Results 
for CP-averaged branching fractions are summarized in table HI table Mil and table [VI while predictions on direct 
CP asymmetries are given in table HI1 table ITVl and table PVTl In B° / B° — > 7r ± / 9 T decays, it is easy to identify 
the final state mesons. Thus one can sum B°/B° — > 7r~ p + up as one channel, although the summed channels 
are not CP conjugates. The B° / B° — > ir + p~ can be summed as another channel and it is also similar for the 
branching ratios of B° / B° — > K*K and B Q S /B® — > K*K decays. In table U and table [V] we give our predictions 
on the summed branching ratios in B°/B Q — > Tr ± p T , K*°K°(K*°K °) and two B s — ► K*K decays. We also give the 
predictions on the sum of the CP-averaged branching ratios of B° — > tt~ p + and B° — > ir + p~ and the other three 
Br a ) decays in table [p and table IVI1 In order to c omp are with QCDF approach 
approach 



46 



47 



48 



42 



43 



44 



49 



50 



51 



53 



54 



3 



451 and PQCD 

55|,l56|,l57|,l58L we also collect their results in these tables, together with 



the experimental data available at HFAG 

Due to several approximations arc made in this work, there are some important possible corrections which we would 
like to address. First of all, our results for the 16 inputs are obtained through the exact flavor SU(3) symmetry for 
the form factors and charming penguins. The amplitudes may receive sizable corrections from the SU(3) symmetry 
breaking effect proportional to m 3 j Kqcd ~ 0.3. Secondly, since we have concentrated on the leading order analysis, 



the radiative corrections proportional to a s (y / mi,AQCD)/"x ~ 0.1 are also neglected. Although we have included one of 
the most important power corrections (chiraly enhanced penguins), the other parts of power corrections proportional 
to A = ^/Agc/j/jTih ~ 0.3 are not incorporated in our analysis. At last, there are also uncertainties from the input 
parameters such as the b quark mass, Wilson coefficients, etc. To characterize these effects, we vary the magnitudes of 
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the non-perturbative charming penguins by 20% and the phases by 20° . We also assume that the gluonic form factors 
C s and £jg have additional uncertainties (±0.05). In the predictions for branching fractions and CP asymmetries 
collected in tables HI Mil IVl ITT1 IIVI and IVI( the first kinds of uncertainties are from these hadronic uncertainties: 
charming penguins and gluonic form factors; the second kinds of uncertainties are from those in the CKM matrix 
elements. 

B. b — > d transitions without 7](') 

b — > d transitions are induced by the operators whose CKM matrix elements are V ub V* d (i = u,c,t). To make it 
clear, we decompose the decay amplitudes into three terms according to the CKM matrix elements: 

A(B - M X M 2 ) = ^|m| {V ub v; d A u + V cb V c * d A c - V tb V t * d A t } , (56) 

where A c is from the charming penguin term. The decomposition is over complete since the unitarity property of 
CKM matrix can be used to eliminate one of the three combinations of CKM matrix elements. We keep all of them 
according to the different dynamics in the corresponding amplitudes. The values for CKM matrix elements: 

\V u bV* d \ = 3.48 x 10- 3 , \V cb V c * d \ = 9.17 x lO" 3 , \V tb V* d \ = 8.60 x 10~ 3 (57) 

will definitely character the branching fractions and CP asymmetries. 

B° — > n ± p If: are dominated by tree operators which has the CKM matrix elements: V ub V* d . To illustrate the 
situation, we will use the second kind of inputs given in Eq. (|54[) and take B° — > p + ir~ as an example (in units of 
GeV): 

\A U (B° -> p + n-)\ = 0.131 x (1.03C y + 0.77CJ ) ~ 260 x 10" 4 , 
\A C (B° -> p+7r-)| = \A™\ ~ (30 ~ 40) x lO" 4 , 

\A t (B° p + n-)\ = |0.131(-0.0015C y - 0.007CJ )| - 5 x 10~ 4 . (58) 

Our predictions on branching fractions of B° — > 7r ± J o =F decays are smaller than those in QCDF [l(3]. Neglecting 
the small terms, the main reason is our smaller B — > P and B — > V form factors: QCDF uses much larger form 
factors F B ^ = 0.28 ± 0.05 and A^ p = 0.37 ± 0.06. In the present framework, BTZ(B Q -> p+vr") is smaller than 
B7Z(B° — ► p~ir + ). In the first solution, the fitted B — > V form factor Aq = 0.233 is almost equal with the B — > P 
form factor F = 0.206. Since the decay constant of p meson is much larger than that of tt: 0.209/0.131 ~ 1.5, we 
expect B7Z(B° — > p + tt~) is only one half of B1Z(B° — > p~ir + ). Charming penguins AY C P and A^Y can slightly change 
the ratio: the charming penguin A^ in B° — > p + ir~ gives a destructive contribution, while A^J* in B° — > p~ir + 
gives a constructive contribution. In the second solution, contributions proportional to form factors are almost equal 
with each other, as the B — ► V form factor Aq^ v = 0.291 is much larger than F B ^ P = 0.198 which can compensate 
differences caused by decay constants. But unlike in the first solution, the role of charming penguin totally changes: the 
charming penguin in B° — > p + ir~ gives a constructive contribution, while A^ p in B° — > p~n + can give a destructive 
contribution. It is reasonable, since the charming penguins A^f and A^Y almost interchanges the phases. 

Our predictions for branching ratios of B° ir° p° are larger than that in QCDF especially the prediction utilizing 
the inputs given in Eq. (|50p . In this channel, two kinds of charming penguin almost cancel with each other, since they 
have similar magnitudes and but different signs as given in Eq. (|50|) and Eq. ([54]) . The tree contribution proportional 
to the soft form factor £ is color-suppressed (the Wilson coefficient Ci + jt- ~ 0.12 is small compared with that of 
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TABLE I: Branching ratios (in units of 10 -6 ) of B — ^ VP decays induced by the b — > d (AS = 0) transition: the first 
solution (This work 1) and the second solution (This work 2). In both cases, we have included the chiraly enhanced penguin 
in B — » VP decay amplitudes. The first kinds of uncertainties are from uncertainties in charming penguins and gluonic form 
factors as discussed in the text; the second kinds of uncertainties are from those in the CKM matrix elements. We also cite the 
experimental data and theoretical results given in QCDF and PQCD 



0, 0, y , 

approach to make a comparison. 



Channel 


Exp. 


QCDF 


PQCD 


This work 1 


This work 2 


B~ -> p"7T° 


10.9+11 

— 1.5 


14 n+e-s+s. 1+1.0+0.8 

— 5.5 — 4.3—0.6 — 0.7 


6-9 


Q Q + 0.3 + 1.0 

— 0.1 — 1.0 


n 4 +o.6+i.i 

— 0.6—0.9 


ET -» pV" 


8.7+1-? 

— 1.1 


-1 -1 n+6. 3+3. 6+2. 5+1. 3 
—5.0—3.1—1.2 — 1.1 


10.4+^i ±2.1 

— 3.4 


ir, y+0. 7+1.0 

— 0.7—0.9 


j q+0. 2+0.8 
— 0. 1 — 0.8 


B~ — ► LU7T~ 


6.9 ± 0.5 


Q q+4.4+2. 6+1. 8+0.8 
— 3.5—2.2—0.9—0.9 


11.3±|-| ± 1.4 


fi 7 +0.4+0.7 
— 0.3—0.6 


o r+0.3+0.8 
— 0.3 — 0.8 


B - _^ k*°K~ 


< 1.1 


n o n +0. 11+0. 12+0. 09+0. 57 

VJ ■ O w A AA A 1 A A AA n 1 rt 

uvv — u.uy — u.iu— u.uy — u.iy 


0-31 + So8 


n ^q+0.26+0.09 
— U.zU — O.Uo 


n ri +0.18+0.07 

U.U 1 A1C A Al? 


B~ -> K*~K° 




n qq+0. 08+0.41+0. 08+0. 58 
U, '- )U — 0.07— 0.18 — 0.07— 0. 17 


i 00+0.68 
—0.47 


a r 4 +0.26+0.10 
— 0.21—0.08 


q ri +0.21+0.08 
—0.17—0.07 


B~ -> (f>TT~ 


< 0.24 


« 0.005 




w 0.003 


w 0.003 


B 1 ' -> p~n+ 1 

B°^p+TV- J 


24.0 ± 2.5 


Of; r + lo.2 + lU.3+2.U+3.y 
" D,0 -14.7- 8.6-3.5-2.9 


18-45 


i q a +0.6+1.2 
i,3 - 4 -0.5-1.2 


1 ^ q+0.5+1.6 
ib - 8 -0.5-1.5 


D j D — > p 7T 






9/1 Q/1 


i 9 a+1.9+1.2 


1/1 Q + l.6+1.5 

14.8.! 5 _ 14 


B / B — » p TT T 






24-34 


1 /l A+1-9+1.3 

14.9I 19 I 13 


1 o t+1.5+1.7 
J-O-l-1.6-1.6 


D — > p TV 




1 c 4+8. 0+5. 5+0. 7+1. 9 
J - J ' fl — 6.4— 4.7-1.3-1.3 




c q+0.5+0.5 
J - J -0.5-0.5 


a R+0.2+0.7 
u - u -0. 1-0.7 


SO — +a 
B — > p 7T 


ion i c\ >-t 

13.9 ± 2.7 


91 o+10. 3+8. 7+1. 3+2.0 
*J-«_ 8.4-7.2-2.3-1.6 




7 r+0.3+0.8 
'•°-0.1-0.8 


1 A O+0.4+0.9 

10.2I 05 I .9 


nO . „0„0 
JD — > p 7T 




a +0.2+0.2+0.9+0.5 
"•^-0.2-0.1-0.3-0.3 


a fi7 a 1 1 
0.U/-U.11 


O r+0.2+0.2 
z ' d -0.1-0.2 


1 r+0. 1 + 0.1 
1 - J -0. 1-0.1 


rO t , ,„0 


< 1.2 


n m +0.00+0.02+0.02+0.03 
u ' ul -0. 00-0. 00-0. 00-0. 00 


A 1 A A OO 

U.1U-U.28 


n nnm+ -o 299 + .oooo 

u.uuuo_ .0000-0.0000 


n m c+o. 024+o. 002 

U.UlO_ ooo-o .002 


rO . 7^*0 S-0 




H oc+0. 08+0. 10+0. 08+0. 46 
u ' ZD -0. 07-0. 09-0. 08-0. 15 




a ^r+0.24+0.09 
u -^°-0.19-0.07 


a ^7+0.17+0.06 
u - 4: ' -0.14-0.05 


r> — > K i\ 


< 1.9 


a 9Q+0. 10+0. 39+0. 08+0. 60 
u ' zy -0.09-0. 17-0.07-0. 17 




a ci +0.24+0.09 
u ' o± -0.20-0.08 


n 4 o+0.20+0.07 
u - ,o -0. 16-0.06 


B° -► iT !? i 






« 1.96 


1 A O A 1 A 1 O 

A Qf;+0.34+0.18 
U.ilO_ .27-0.15 


1 A 1 A 1 A 

A AC+0.26+0.14 
U.»0_o.22-O.12 


B°/B°^K*°K 








A Qir+0.34+0.18 
u ' 3J -0.27-0.15 


Q4+0.26+0.14 
u ' i,y: -0.22-0.12 


B /B°^K*°K 








n q 7 +0.35+0.18 
U.»( -0.27-0. 15 


n Q7 +0.26+0.14 
u - a ' -0.22-0.12 


rO j, 

i> — > 07r 


< 0.28 


« 0.002 




W 0.001 


« 0.001 




5.4 ±1.2 


Q ,+4.6+3.6+0.7+0.7 
y -^-3.7-3.0-0.4-0.7 


o r+3. 0+0. 8+0. 4+1. 2 6 
°-°-2. 1-0.7-0.4-0.2 


q Q+2.0+0.4 

°- 3 -1.7-0.4 


q q + 1.9 + 0.3 

•J-* 5 -!. 6-0.3 




q 1+3.7 


R q+3. 1+2. 4+0. 5+0. 5 
u -°-2.5-2.0-0.3-0.5 


o y+3. 0+0. 7+0. 5+1.1 6 
°- ' -2.2-0.9-0.7-0.3 


a q 7 +2.46+0.08 
u - ' -0.22-0.07 


n ^^+3.18+0.06 
'-'■^-0. 20-0. 05 


73° ^ p ^ 


< 1.5 


a n o+0. 02+0. 16+0. 02+0. 05 
U.UO_g .01-0 .10-0 .01-0 .02 


n AO/l +0.012+0. 004+0. 002+0. 102 6 
"■ uz *-0. 007-0. 002-0. 002-0. 005 


a p. A +0.20+0.00 
u - uq: -0.01-0.00 


n 1^+0.33+0.01 
u - iy: -0.13-0.01 


s° -» pV 


< 1.3 


a m +0.01+0.11+0.02+0.03 
U.U1_ oo-0 .06-0 .00-0 .01 


fl n«1 +0.030+0.004+0.003+0.114 6 
U.UOl_ .018-0.003-0.003-0.009 


a /iq+2.51+0.05 
U.4O_ 12 _o .05 


1 n +3. 5+0.1 
1M -0. 9-0.1 


73° -> lj?7 


< 1.9 


n oi +0.14+0.16+0.35+0.22 
U - OJ --0. 12-0.11-0.14-0. 16 


u - z '-0.10 


n Q1 +0.66+0.09 
u - yl -0.49-0.09 


1 ,+0.8+0.1 
J-^-0. 6-0.1 


73° UJT)' 


< 2.2 


n on+0- 10+0. 15+0. 25+0. 15 
U.ZU_ .08-0.05-0.10-0.11 


U.U/O_ .033 


a lo + l. 31+0.04 
u - lo -0.10-0.03 


q 1 +4.9+0.3 
°- i -2.6-0.3 


B° -> 077 


< 0.6 


« 0.001 


63+ 0033 

U.UUOO_ .ooi9 


w 0.0004 


w 0.0008 


73° -» 0??' 


< 0.5 


« 0.001 


73+ 0035 

u - uu ' °-0.0026 


» 0.0001 


w 0.0007 



a We quote the branching ratios for 73 — + p+7r and B ' — > p+7r from_Ref. |59| . 

'For B —> pr) decays, there are two different predictions given in Ref. 15 211 according to the different mixing angles between rj and rj ' . We 
quote the results in which 9p = —10° is used. There are not too many changes for the other predictions as the value for the mixing angle 
dp = —17° is very close to the first one. 



B° — > p+ttT Ci + j^- ~ 1.03), thus the branching fractions of B° — > 7r°p° in QCDF approach and PQCD approach 
are much smaller than B1Z(B° — > p ± 7r T ). One important feature of the SCET framework is: the hard-scattering 
form factor Q is relatively large and comparable with the soft form factor £. Besides, this term has a large Wilson 
coefficient b{, since C2 ± -^-(1 — tjOCi ~ 1-23 is large, it can give larger production rates which are consistent with 
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TABLE II: Direct CP asymmetries involving 6 — > d (AS = 0) transitions: the first solution (This work 1) and the second 
solution (This work 2). In both solutions, we have included the chiraly enhanced penguin in B — > VP decay amplitudes. The 
first kinds of uncertainties are from uncertainties in charming penguins and gluonic form factors which are discussed in the 
text; the second kinds of uncertainties are from those in the CKM matrix elements. We also cite the experimental data and 
theoretical results given in QCDF [nj and PQCD Q, 0, I54I [5?! approach to make a comparison. 



Channel 


Exp. 


QCDF 


PQCD 


This work 1 


This work 2 


B~ -> p-jv 

B~ -> P °TT- 

B~ — > ujn~ 
B - _> k*°K- 
B~ -> K*~K° 


2 ±11 
7+12 

'-13 

-4 ±6 


, n +l. 2+1. 8+0. 4+17. 5 
^• u - 1.2-2. 2-0. 4- 17. 7 
, 1+1.3+2.2+0.6+19.0 
^■-"--0. 9-2. 0-0. 7-18. 8 
1 n+0.5+2.7+0.8+2.1 
i -°-0.5-3.3-0.7-2.2 
9 o r+6. 9+7.8+5. 5+25. 2 
zo -°-5.7-9.0-6.5-36.8 
10 ,,+3.7+7.8+4.2+27.4 
10i *- 3. 0-3. 5-4. 7-36. 7 


0-20 
-20-0 
~ 
-20±5±2 


1 c r + 16.9+1.6 
— 18.9—1.4 

If! S+13.1+0.9 
lu -°-12 .7-0 .7 

n r + 19. 1 + 0.1 
u - d -19. 6-0.0 
q R +6.1+0.4 
°- u -5.3-0.4 

i c+2. 6+0.1 
L -°-2. 3-0.1 


19 q+9.4+0.9 
lz -°-10. 0-1.1 

iq 9+15.5+1.7 
IJ,Z -13.4-1.9 
9 q+13.4+0.2 
^•°-13. 2-0.2 
_a ,+4.1+0.2 
^•^-4.1-0.2 
1 9+I.7+O.I 
x -*—l. 7-0.1 


B° p+n- 
B° p-n+ 
B° p%° 
B° a;7r° 
B° -» K*°K° 
B° f*°]f 


-18 ±12 
11 ±6 

-30 ± 38 


n c+0. 2+1. 3+0. 1+11. 5 
u -°-0.1-1.6-0.1-11.7 
a r+0. 4+1. 2+0. 2+8. 5 
J-- -!). 4-1. 3-0. 3-8. 4 

1 c 7+4.8+12.3+11.0+19.8 
10- < -4.7-14.0-12.9-25.8 

9fi 7+7.4+7.2+5.7+10.9 
zu - 1 -5.7-9.0-6.9-13.4 
a n , +3.8+5.4+4.5+5.8 
-3.0-2.9-5.2-7.4 


-75-0 
-20-75 


q q+17.2+0.9 
J - J -16.7-0.7 
i i + 17.5+1.2 
1J --°-20. 0-1.1 

n fi +21. 4+0.1 
U - D -21. 9-0.1 
_Q 4 +24.0+l.l 
a ' a -0.0-0.9 
q c+6. 1+0.4 
°- u -5.3-0.4 
, r+2. 6+0.1 
L -°-2. 3-0.1 


19 ,+17.6+1.1 
lz -*-15.3-1.2 

in O+9.4+0.9 
J.U.O_ 10 .2-1.0 

q r + 21.4+0.3 
°- d -20.3-0.3 
q Q r+79.1 + 3.4 
«J». J _185. 5-3.1 
_ 4 1+4.1+0.2 
^•^-4.1-0.2 
1 9+I.7+O.I 
J-^-l. 7-0.1 


B~ — > p~r; 


1 ± 16 
-4 ±28 


9 ,+0.7+6.3+0.4+0.2 
^■^-0. 7-6. 3-0. 4-0. 2 
, 1+1.2+7.9+0.5+7.0 
*' A -1. 1-6. 9-0.8-7.0 


1 Q+l.2 + 2 

±J -0.5-14 

i O+3.0+1 
- 18 -1.6-14 


fi fi +21.5+0.6 
u - u -21.3-0.7 
in 0+66.5+2.8 
1J -°-37. 5-3.1 


Q 1+16.7+0.9 
a - x -15.8-0.8 
91 7+135.9+2.1 
- 1 - ' -24.3-1.7 


B° A 
B° - A' 

S° -f wij' 




no /I + 10. 0+65. 3+20. 9+19. 2 
OJ - - 9.5-55.8-21.4-20.8 
n 9+0.1+53.0+11.6+19.2 
u - z -0. 1-76. 5-11. 5-20.1 


1 Q+l.2 + 2 

±J -0.5-14 

1O+3.0+1 
~ i8 -1.6-14 

-69.i±i§:i 
13.9H:; 


A(\ 7+170.4+2.9 
w - 1 -74.3-3.7 
ri 7+103.3+3.4 
J±i ' -42.9-3.9 
q 1+30.7+0.9 
"•^-30.2-1.0 
AO n +87.5+4.8 
*°- u -38 .8-5 .1 


qq q+66. 9+3.1 
JO -°-62.4-2.8 
r 9 9+19.9+4.4 
J ^- z -80. 6-4.1 
q c+17.8+0.9 
y - u -16.8-0.9 
97 9+I8.I+2.4 
z ' - z -29.7-2.2 



the present experimental data. The agreement is very encouraging. 

Branching ratios of B — ► A'* A are larger than those in QCDF for the presence of charming penguins. In B~ — > 
A* - A' and S — > A* AT , both of penguin operators and charming penguins can give contributions. The difference 
for these two channels is: the spectator antiquark in B~ — ► K*~K° is u and it is d in B° — > A* A . It does not 
affect the contributions from either penguin operators or charming penguins, thus we expect the relations B1Z(B~ — > 
A*-A°) = BK(B° -> A* A ) and ^ CP (S^ -> A*"A°) = A CP ( J B° A* A ). The small differences in branching 
fractions are induced by the different lifetimes of _B~ and -B . The analysis is similar for the other two b — > modes: 
S" -> A'- A* and 5° -> A°A*°. 

For the decays with sizable branching fractions, our predictions on direct CP asymmetries are typically small and 
most of them have the correct sign with experimental data. Predictions in QCDF approach on these channels are 
also small in magnitude, but some of them have different signs with our results and experimental data. In PQCD 
approach, the strong phases mainly come from the (S — P)(S + P) annihilation operators. These operators are chiraly 
enhanced and the imaginary part are dominant. Thus the direct CP asymmetries in PQCD approach are typically 
large in magnitude. 
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TABLE III: Branching ratios (in units of 1(T 6 ) for As = 1 processes: the first solution (This work 1) and the second solution 
(This work 2). In both solutions, we have included the chiraly enhanced penguin in B — > VP decay amplitudes. The first 
kinds of uncertainties are from uncertainties in charming penguins and gluonic form factors which are discussed in the text; the 
second kinds of uncertainties are from those in the CKM matrix elements. We also cite the experimental data and theoretical 
results given in QCDF [l(J and PQCD [51L to make a comparison. 



Channel 


Exp. 


QCDF 


PQCD 


This work 1 


This work 2 


B~ 


-> K*~w° 


6.9 ±2.3 


0+1. 1+1.0+0.6+4.4 

"•°-1.0-0.9-0.6-1.4 


4-3i 5 2 ;° 2 


/l 9+2.2+0.8 
^•^-1.7-0.7 


e c + 1. 9+0. 7 
u -°-1.7-0.7 


B~ 


K*°n- 


10.7 ±0.8 


q R +0. 4+1. 5+1. 2+7. 7 
0.3-1.4— 1.2-2.3 


6.ol?;? 


o (.+4.7+1.7 
°- J -3.6-1.4 


Q q+3.5+1.3 
y ' a -3. 0-1.1 


B~ 


-> P °K- 


4.25tS;g 


9 e+0. 9+3. 1+0.8+4.3 
^• u -0. 9-1. 4-0. 6-1. 2 


5-1-2.8 


„ +2.7+1.0 
u - ' -2.2-0.9 


4 fi +l. 8+0.7 
4 ' D -1.5-0.6 


B~ 


-^p-R° 


s.oiii 


r o+0. 6+7. 0+1. 5+10. 3 
J '°-0.6-3.3-1.3- 3.2 


8.7_ 4 ; 4 


+4.7+1.7 
y -°-3.7-1.4 


lf) ,+4.0+1.5 
lu -J-3.3-1.3 


B~ 




6.7 ±0.5 


o c+1. 0+3. 3+1.4+4.7 
"•"-1.0-1.6-0.9-1.6 


10-6_ 5 ° 8 4 


c 1+2.4+0.9 
"• J--1. 9-0.8 


r q+2. 1+0.8 
"• a -1.7-0.7 


B~ 


-> <j>K~ 


8.30 ±0.65 


a r+0. 5+1. 8+1. 9+11. 8 
^•"-0. 4-1. 7-2.1- 3.3 




q 7+4.9+1.8 
J - 1 -3.9-1.5 


Q R +3.2+1.2 
°' u -2. 7-1.0 


B° 




o.o±i;? 


n 7+0.1+0.5+0.3+2.6 
u - ' -0.1-0.4-0.3-0.5 


2-0_ j 


a R +2.3+0.9 
*-°-l. 8-0.7 


q 7+1.4+0.5 
°- ' -1.2-0.5 


B° - 


-> K*-TT + 


9.8 ± 1.1 


q q + 1. 4+1. 3 + 0. 8 + 6. 2 
"•"—1.1-1.2-0.8—1.6 




o ^+4.4+1.6 
8-4-3.4-1.3 


Q r+3.2+1.2 
V-°-2. 8-1.1 


B° 


-» p°K° 




a K +0. 5+4. 0+0. 7+6.1 
^• u -0. 5-2. 1-0. 7-2.1 




q c+2.0+0.7 
°- J -1.5-0.6 


r O+2.1+0.8 
J -°-1.8-0.7 


B° 


-> p+R- 




7 4 + 1.8+7.1+1.2+10.7 
'•^-1.9-3.6-1.1- 3.5 


8.8tl s 5 


O Q+4.6+1.7 
3 -°-3.7-1.4 


i n 9+3.8+1.5 
±u - z -3.2-1.2 


B° 




5.0 ±0.6 


9 q+0. 3+2. 8+1. 3+4. 3 
^•°-0. 3-1. 3-0.8-1. 3 


9.8±H 


a t +2.1+0.8 
^•J-l.7-0.7 


4 q+1. 9+0.7 
^• a -1.6-0.6 


B° 




8 3+ 1 ' 2 
°-°-i.o 


a 1+0.4+1.7+1.8+10.6 
^^-0. 4-1. 6-1. 9- 3.0 


7-3±f;l 


n 1+4.6+1.7 
y -J-3.6-1.4 


o q+3. 0+1.1 
°' u -2. 5-1.0 


B~ 


-> K*-r] 


19.3 ± 1.6 


111 Q+l.9+8.1+1.8+16.6 
J-U.o_l.7_4.4_l.3_ 5.5 


99 i q+0. 26 
__.lO_ 27 


1 7 q+5.5+3.5 
1 ' - y -5.4-2.9 


IS R+4.S+2.5 
JO.O_4 g_2 2 


B~ 


-> K*~T)' 


4 q+2.1 


r i +0.9+7.5+2.1+6.7 
0-3.8-3.0-3.3 


6.38 ± 0.26 


a r+6.6+0.9 

^•^-s.g-o.s 


a cj+5.3+0.8 
4-8_3. 7 _o.6 


B° 


-» K*°T] 


15.9 ± 1.0 


Ifi 7+1-1+7.8+1.4+16.2 
±Ui ' -1.0-4.3-1.2- 5.5 


99 qi +0.28 
__.ol_ 29 


lfi fi +5.1+3.2 
io -°-5.0-2.7 


16 S+ 4 ' 1+2 ' 3 
10 - J -4. 3-2.0 


B° 


A*°77' 


3.8 ± 1.2 


q q+0. 4+6. 6+1. 8+6. 2 
°- a -0.4-3.3-2.5-2.9 


q qr + 0.29 

°- JJ -0 .27 


a 1+6.2+0.9 
*- A -3. 6-0.7 


4 n +4.7+0.7 
4 - u -3.4-0.6 



C. b — > s transitions without r; and 77' 

Like b d processes, 6 — > s decay amplitudes can also be decomposed into three different parts according to the 
CKM matrix elements. The values of the CKM matrix elements are given by: 

\V ub V* s \ = 0.81 x 10- 3 , \V cb V*\ -39.41 x LCr 3 , |tfcV£| = 40.66 x 10" 3 . (59) 

Tree operators are highly CKM-suppressed, but the CKM matrix elements for the rest two kinds of contributions A c 
and A t are in similar size. Together with the hierarchy in Wilson coefficients: „> C^-iO) charming penguins will 
provide a dominant contribution. For example, the penguin operators in B~ — ► tt~K° decay process is proportional 
to (24 + r x ae, B~ — > Tr^A"* is proportional to 04 while B~ — ► p~K° is proportional to (14 — r x ag, where 04,6 = 
(74,6 + Cz$/N c and r x = 2/xp/m;,. Thus if we only consider the emission diagrams, B1Z(B~ — > 7r~.A ) > B7Z(B~ — > 
7r~A'*°) > B1Z{B~ — > p~ K°) holds, since 04 ~ a6 and r x ~ 1. But in the present framework, contributions from 
penguin operators proportional to VtbV t * do not play the most important role: 



1^(5- 




= [0.16 x (-0.044C P - 0.036Cj )| - 


15 x 10~ 4 , 






= |0.16 x (0.0004C y ± 0.004CJ )| ~ 


1 x 10~ 4 , 






| = |0.217 x (-0.022C P - 0.015Cj )| 


- 10 x 10" 4 



Compared with the results given in Eq. ([50]) and Eq. (f54l) . we find penguin operators are smaller than charming 
penguins. According to the size of charming penguins, we expect the relation B1Z(B~ — > p~K°) ~ B1Z(B~ — > 7r _ A'* ). 
This is well consistent with the experimental data. 
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TABLE IV: Direct CP asymmetries (in %) for As = 1 processes: the first solution (This work 1) and the second solution (This 
work 2). In both solutions, we have included the chiraly enhanced penguin in B — > VP decay amplitudes. The first kinds of 
uncertainties are from uncertainties in charming penguins and gluonic form factors which are discussed in the text; the second 
kinds of uncertainties are from those in the CKM matrix elements. We also cite the experimental data and theoretical results 
given in QCDF ljj and PQCD [51L to make a comparison. 



Channel 


Exp. 


QCDF 


PQCD 


This work 1 


This work 2 


B~ -> K*~tv 


4 ±29 


7+2.1+5.0+2.9+41.7 
°- ' -2.6-4.3-3.4-44.2 


09+21 
°^-28 


1 7 O+30.3+2.2 
1 ' •"-24.6-2.0 


i 9 q+12.0+0.8 
u ' a -12.2-0.8 


B~ -> K*°iv- 


-8.5 ±5.7 


1 fi+0. 4+0. 6+0. 5+2. 5 
^•"-0. 5-0. 5-0. 4-1.0 


-l±o 








B~ -> p°K~ 


01+11 

Ji -10 


10 c+4. 5+6. 9+3. 7+62. 7 
J- - - 5. 7-4.4-3. 1-55.4 


7lli 


Q 9+15.2+0.7 
y - z -16.1-0.7 


, R n +20.5+1.3 
10.U_22.4-1.6 


B- -» p-K° 


-12 ± 17 


n o+0. 1+0. 3+0. 2+1. 6 
u -°-0. 1-0.4-0. 1-1. 3 


1 ± 1 








B~ -* cjR- 


2±5 


7 o+2. 6+5. 9+2. 4+39. 8 
'■°-3. 0-3. 6-1. 9-38.0 




i, fi +18.2+l.l 
11 - u -20. 4-1.1 


i 9 O + 16.6+0.8 
— 17.3—1.1 


B~ (f>K~ 


3.4 ±4.4 


i R +0. 4+0. 6+0. 5+3.0 
1 -°-0. 5-0.5-0. 3-1. 2 


1+0 
1 -i 








B° -> K*%° 




19 B+4. 0+4. 7+2. 7+31. 7 
lz -°-3.2-7.0-4.0-35.3 




r n +7.5+0.5 
J - u -8.4-0.5 


r 1+4.8+0.4 
°- 4 -5. 1-0.5 


B° -> A'*-7r+ 


-5 ± 14 


n 1+0.6+8.2+5.1+62.5 
z>i -0.7-7. 9-5. 8-64. 2 


-eotg 


i i 9+19.0+1.3 
11-^-16.2-1.3 


i 9 9+II.4+O.8 
iz - z -11.3-0.8 


B° -> p°K° 


-2±27±8±6 


7 r + 1. 7+2. 3+0. 7+8. 8 
' -°-2. 1-2. 0-0. 4-8. 7 


7+s 
'-5 


R R +11.6+0.8 
u - u -9. 7-0.9 


n r + 4.8+0.3 
J -°-4.8-0.2 


B° -» p+A" 


22 ± 23 


o o + l. 3+4. 4+1. 9+34. 5 
J '°-1.4-2.7-1.6-32.7 




7 i+H.2+0.7 
' • 1 -12.4-0.7 


Q c + 13.0+0.7 
a - u -13.5-0.9 


B° -> c;A 


21 ± 19 


o 1 +2.5+3.0+1.7+11.8 
°' 1 -2.0-3.3-1.4-12.9 


0+2 
J -3 


r 9+8.0+0.6 
J - z -9.2-0.6 


n q+5.2+0.3 
°-°-5. 4-0.3 


S° -» ^>A° 


1 ± 12 


i 7+0.4+0.6+0.5+1.4 
L -' -0.5-0.5-0.3-0.8 


o+l 

°-2 








B" -> A*"?? 


2±6 


o r+0. 9+1. 9+0. 8+20. 7 
J - d -0.9-2.7-0.8-20.5 


-24.57lS;g 


9 R +5.4+0.3 
^•"-5.5-0.3 


i q+3. 4+0.1 
x - a -3. 6-0.1 


S" -> K*~rj' 


30±P 


■> A 9 +4.7+ 8.5+ 4.9+27.5 
14: - z -4.2-13.8-14.6-26.1 


4.60^ 


9 7+27.4+0.4 
Z -'-19. 5-0.3 


9 R +26.7+0.2 
z - D -32.9-0.2 


B° -» A*°r7 


19 ±5 


o q+0. 9+1. 1+0.2+3. 8 
— 1 . 1 — 0.8 — 0.2 — 3.5 


0.57 ±0.011 


i 1+2.3+0.1 
1 - i -2. 4-0.1 


n 7+I.2+O.I 
u - ' -1.3-0.0 


B° -» K*°rf 


-8 ±25 


r r + 1. 6+3. 1 + 1. 8 + 6. 2 
°-' J -l. 3-5. 1-5. 9-7.0 


-1.30 ±0.08 


q fi +8.9+1.3 
"•"-11.0-1.2 


Q q+6.2+0.9 
y - y -4.3-0.9 



From table HVl we can see the direct CP asymmetries of B~ — ► K*°tt~ , B~ — > K°p~ , B~ — ► K~4> and B~ — ► A"°(/> 
are zero. In these channels, tree operators do not contribute. The weak phases for penguin operators and charming 
penguins are equal to each other, which can not induce any direct CP violations. CP asymmetries in other channels 
are not large, because the strong phases of charming penguins are either close to 0° or 180° and imaginary parts are 
accordingly small. The PQCD results for most B — > K*ir and B — ► pK channels are much larger than ours, since 
they have more large imaginary part from annihilation diagrams. The QCDF results are small and comparable with 
ours but with a relative minus sign. We have to wait for the experiment data to resolve this disagreements. 



D. B Decays involving r\ or rf 

As we can see from table U there is about 3.1a deviation for our prediction on the branching ratio of B~ — > p~r}' 
from the experimental data. Contributions from penguin operators are suppressed by the CKM matrix elements 
as given in Eq. (|57|) and the dominant contribution is from the tree operator. This kind of contribution is either 
proportional to B — > rj q or B —> tj s form factor. Utilizing results given in Eq. (f5U|) and Eq. (|54"|) . we obtain B — > r\ q 
and B — > r) s form factors as follows: 

pB^ nq = fcP + + 2 Q + 20 a ) = (0.053 ± 0.068) [(0.100 ± 0.021)] , 

pB—>7} s = (Q + (; Jg ) = (-0.076 ±0.055) [(-0.049 ± 0.011)], (61) 
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where the results in (out) the square brackets are predictions using the second (first) kind of inputs. In equation (|61[) . 
we can see: after taking the gluonic form factors into account, the F B ^ r)q and F B ^ Vs form factors are in the similar 
size but with different signs in both kinds of inputs. In B~ — > p~ r) q , another tree operator contributes in which rj q is 
emitted. Although this contribution is color-suppressed, terms proportional to (J give a sizable contribution. It can 
be estimated by using a larger effective B — > r) q form factor. Recalling that physical states 77 and rj are mixtures of 
rj q and rj s as in Eq. (|22[) . one obtains the expressions for B — > rj(') form factors: 

v2 

A s ^' = sin(<9) + A B ^" 3 cos(0). (62) 

v2 

The mixing angle between ?7 g and ry s has been determined as 9 = (39.3 ± 1.0)° 33j, [3J, [35| which is very close to 45°, 
thus we can obtain very small B — ► ?/' form factors and relatively large B — > rj form factors. Thus the branching 
fraction of B~ — > p~if is relatively suppressed for this flavor structure. In QCDF and PQCD approaches, the form 
factors are different: F B ^ Vq ^ F ne _ Thus the predicted branching ratio of B~ — > is comparable with 
B1Z(B~ — > p^ry') in these two approaches. 

As in B° — > 7r°p° process, our predictions on branching fractions of B° — > p°rj(') and _B° — > u^C) are much larger 
than the results evaluated in QCDF and PQCD approach. These channels are the so-called color-suppressed decays, 
as the contributions from terms proportional to £ and Q are small due to the small Wilson coefficients. But in the 
present framework, the hard-spectating form factors (j and £j ff are comparable with £ and £ g . Moreover, the Wilson 
coefficients for these form factors are large. Thus branching ratios of B° — > p°?y(') and B° — > lutjC) are much larger. 

Similar with B — > K*tt and B — ► pK decays, B — > K*tj(t]') are also induced by 6 — ► s transitions in which charming 
penguins provide most important contributions. But compared with £> — > A^*7r and _B — > pA decays, there are 
something new in these channels. In B — > K*r){rj') } there exist three kinds of charming penguins: 

^cc Vq = 7/1^"^ ^ccs)' 1,3 = + A^, v . (63) 
Substituting the values given in Eq. (|50p and Eq. (|54p . we obtain ratios of charming penguins: 

|^jft(A^ + 2Ag) - ringKjg + Af/)| 2 2 q 



(A^ + + sin(0)(A^ + 

The branching fraction of B° — > A*°?7 is about 4 times larger than that of B° — > K*°r]' for both solutions. The main 
reason for the difference is: A^ 1,3 is very small due to the cancelations between A^ y and A^; the penguin operators 
play the dominant role in the B — > K*r\ s decay amplitudes. Our results for these channels have a better agreement 
with experiments than QCDF and PQCD. 



E. B s VP Decays 

Since we have assumed the SU(3) symmetry for form factors and charming penguins, branching fractions and direct 
CP asymmetries of the B s decays are related to the corresponding B decays: 

B1Z(B° S -> K* + K~) = BTZ(B° -» p + K~), BTZ{B° S -> K + K*~) = BTZ(B a -> tt+A*"), (64) 
A CP (B° S -> A*+A-) = A C p(B° -> p+A-), A CP (B S ° -> A+A*~) = A C p(B° -> tt+A*-). (65) 
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These relations can also be applied to the following channels: 

BK{B° S -> A'* + tt") = BTZ{B° -> p+tt"), BK(B° S ^ K + p~) = BTZ(B° -> 7T+/T), (66) 
A C p(5 s ° ^ A*+0 =^cpOB° ^ P+O, A CP (B° - tf+p") = AcpCS -h. 7T+p-). (67) 

In tree-operator-dominated processes _B° — > p~ K + , we obtain branching ratios which are much smaller than 
predictions in the other two approaches: because PQCD predicts F Bs ^ K = 0.24^ ; 04 ^ ; 01 and QCDF use an even 
larger form factor F B ^ K = 0.31 ±0.05. BTl(B° -> ir~ K*+) is consistent with results in QCDF and PQCD approach 
as the B — > K* form factors are consistent. As in B decays, we also predict larger branching ratios for color-suppressed 
B s decays than QCDF and PQCD which can be tested on the future experiments. 

Our predictions onfus processes B® — > K* K are consistent with the other two approaches. But there are huge 
differences in our predictions of B7Z(B S — * (frviv')) with those in QCDF and PQCD. In PQCD approach, contributions 



from gluonic components of 77 and 7/ in B — > t]{') form factors are very small and can be neglected 60(. As shown in 



Ref. 



decay amplitudes of B s — > <pr] q are dynamically enhanced sizably, as the Wilson coefficients 03 — strongly 
depend on the factorization scale. In B s — > 4>r) s , dominant penguin operators are either proportional to 04 — 2r x ae or 
04. The former Wilson coefficient is very small as 04 ~ etg and 2r x ~ 1. The total decay amplitudes of B s — ► <prj q and 
B s — > (f>rj s are in similar size but with different signs. Thus branching ratio of B s — > 4>rj predicted in PQCD approach 



is relatively large while branching ratio of B s — > <prj' is small due to cancelations between the two amplitudes [58| . In 
the SCET framework, charming penguins play the most important role: the charming penguin AY C P almost cancels 
with A^Y ■ Thus the dominant contributions to B s — > (j>rj(rj') are from the gluonic charming penguin and the penguin 
operators which are proportional to VtbV t *. Neglecting the latter term, we have: 

A*r^ = co S (6)V2Al P g sm(8)Al P ~ (V2 - l)A v cc P , 

A P r H ' = sM0)V2A v c p + cos(6)(Al P ) ~ (V2 + l)A v c p . (68) 

These two equations can explain the small branching fraction for B s — > 4>rj together with the large one for B s — > (£>r]' . 
The QCD penguin contributions do not change the ratios too much, but sizable differences appear in the two solutions. 
The large differences in two kinds of predictions on direct CP asymmetries also confirm this feature. 

In B s decays, there are 7 decays in which the direct CP asymmetries are zero: B s — * K*°K°, B s — > K*°K°, 
B s — > 7T°(p and B s — > p a (ui)ri('). As we know, in order to give a non-vanishing direct CP violation, at least two decay 
amplitudes with different weak phases and different strong phases are required. In the first two decays, contributions 
from tree operators vanish at leading order. The non-zero contribution is either proportional to the CKM matrix 
elements VtbV t * s or V c bV* s and both of them are taken real in our calculation. Thus in these two channels, there are 
only one weak phase and direct CP asymmetry is in the present framework. The latter 5 channels are induced by 
b — > s transitions and one of the final state mesons is neither open nor hidden strange. There is no contribution from 
charming penguins in these modes. The direct CP asymmetries are zero for lack of necessary strong phases. 

F. Mixing-induced CP asymmetries 

In this subsection, we will discuss mixing-induced CP asymmetries which can be studied via time-dependent mea- 
surements of decay widths. The four decay amplitudes in B° / B° — ► /(/) decays are defined by: 

A f = (f\H eff \B°),A f = (f\H ef f\B°),A f = (f\H eff \B°), A f = (f\H eff \B°). (69) 
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TABLE V: CP-averaged branching ratios (xlCP 6 ) of B B — > PV decays: the first solution (This work 1) and the second solution 
(This work 2). In both solutions, we have included the chiraly enhanced penguin in B — > VP decay amplitudes. The first 
kinds of uncertainties are from uncertainties in charming penguins and gluonic form factors which are discussed in the text; 
the second kinds of uncertainties are from those in the CKM matrix elements. We also cite theoretical results evaluated in 
QCDF [3 and PQCD [3] to make a comparison. 



Modes 


QCDF 


PQCD 


This work 1 


This work 2 


B° -» K + K*~ 
B° s -> K*+K~ 
B° _» K°K*° 
B s ° -» K*°K° 
B° s /B° - K+K- 

B ° s /B° -» K*+K- 
B° — > -i 

B° a ^K*°K° > 

m -> rrv 


j - 7_|_1 D-J-Q 9 

1. 5—1. 3-0. 9-2.3 
r r + 1. 3+5. 0+0. 8+14. 2 
J - J -1.4-2. 6-0.7- 3.6 

q Q+0. 4+1. 5+1. 3+10. 4 
■J-»_o.4_l, 4-1.4- 2.8 

a 9+0. 4+4. 6+1. 1 + 13. 2 
^■ z -0.4-2.2-0.9- 3.2 

n i n+0.03+0.04+0.01+0.02 
u - iz - 0.02-0. 04-0. 01 -0.01 


„_i_i 7_|_1 7_|_n 7 

°- u -1.5-1.2— 0.3 

A 7+1.1+2.5+0.0 

' -0.8-1.4-0.0 
7 0+2. 5+2. 1+0.0 
' 1.7-1.3— 0.0 
1 q+0. 7+2. 2+0.0 
*-°-0. 7-1. 4-0.0 

n 1 R+0- 06+0. 02+0. 00 
U ' 1D -0. 05-0. 02-0. 00 


O 0+4.6+1.7 
3 '°-3.7-1.4 

7 q+4.4+1.6 
' •<*— 3.4-1.3 

7+4.4+1.6 
°-'-3.5-1.4 

Ifi K+6.4+3.2 
10.0_ 4 g _2 g 

iq g+6.9+3.4 
la -°-5.6-2.9 

1 q 9+6.3+3.3 
io - z -5.0-2.7 

Ifi c+6.2+3.2 
10.0_ 4 9 _ 2 7 

n n7 +0.00+0.01 
u - ul -0.00-0.01 


-J. 1 } 9_|_1 9 

J -°-2. 8-1.1 
in 9+3.8+1.5 
lu - z -3.2-1.2 

q q+3.2+1.2 
a -°-2. 8-1.0 

1+3.7+1.4 
y -*- 3.1-1.2 

1 7 r+5.0+2.5 
1 ( .J_4.4_2.1 

91 0+5.4+2.8 
z J--8_4.7-2.4 

1 q 7+5.O+2.6 
!»■ ' -4.2-2.2 

1 « 7+4.9+2.6 
IO. < _4.2-2.2 

n fiQ+o oo+o o 1 

U.U3_ 00— 0.01 


B° s ^p°K 
B s ° -f 


o 7+4.6+3.5+0.7+0.8 
°- ' -3.7-2.9-1.0-0.7 

n 9C+0. 08+0. 10+0.32+0.30 
U.ZO_ .08-0 .06-0 .14-0 .14 

91 c+11.9+9.2+1.8+1.6 

zy: - d -9.7-7.8-3.0-1.6 
n R1 +0.33+0.21 + 1.06+0.56 
u - Di -0. 26-0. 15-0. 38-0. 36 
n ri +0.20+0.15+0.68+0.40 
u - di -0.18-0. 11-0. 23-0. 25 

n 97+0.09+0.28+0.09+0.67 
u - z ' -0.08-0.14-0.06-0.18 


7 R +2.9+0.4+0.5 
' - u -2.2-0.5-0.3 
n n7 +0. 02+0. 04+0. 01 
u - u '-0. 01-0. 02-0. 01 
17 0+7.7+1.3+1.1 
J -'' o -5.6-1.6-0.9 

n ns+0-0 2 +0-0 7 +0- 01 

U.UO_ . 02-0. 03-0. 00 
n 1 r+0. 05+0. 07+0. 02 
u ' ld -0. 04-0. 03-0. 01 

n 1C+0. 04+0. 09+0. 02 
u - 1D -0. 03-0. 04-0. 01 


r q+0.5+0.5 
°- a -0.5-0.5 

n an+o oT+O- 10 

U - MU -0.01-0.11 
7 R +0.3+0.8 
' '"-0. 1-0.8 
9 q+0.2+0.2 
z - u -0.2-0.2 
n Qn +0.08+0.10 
u - JU -0.01-0.11 

n 11+0.23+0.08 
"■^-0.18-0.07 


fi fi +0.2+0.7 
u - u -0. 1-0.7 

1 Q7+0.16+0.10 

J^-u'-o.is-o.og 
ln 9+0.4+0.9 

iu - z -0.5-0.9 

n 01 +0.05+0.08 
U.Ol_ .02-0.09 

1 q + 0. 1+0.1 

^•"-0. 1-0.1 
n ri+0.21+0.08 
u - dy: -0. 17-0.07 


B° s - A 
B° -> a, ?7 


n 17+0.03+0.07+0.02+0.02 
u - i ' -0.03-0.06-0.02-0.01 
n 9C+0. 06+0. 10+0.02+0.02 
U.ZO_ .05-0.08-0.02-0.02 

n nn+0- 005 + ' 01 o+0'0 28 +0'0 25 

U.U1Z_ 004 _o. 003-0. 006-0. 006 

n no/i+0- 011 +0- 028 +0- 077 +0' 042 

U.UZ4_ 009-0.006-0.010-0.015 

n 19+0.02+0.95+0.54+0.32 

u - iz -0. 02-0. 14-0. 12-0. 13 
n qe+0. 01 + 1. 10+0.18+0.40 
U-UO_ o 01 -0.17-0. 08-0. 04 


n nR+ -o 3 + - 01 + -oo 
u - uu -o. 02-0. 01 -0.00 

n 1 q+0. 06+0. 02+0. 00 

U.lO_ . 04-0. 02-0. 01 
n n i +0.03+0.05+0.00 
u - u4 -0. 01-0. 02-0. 00 
n a 1+0.18+0.15+0.00 
u -^—0. 13-0. 14-0.01 
q R +1.5+0.8+0.0 
•J-O-l. 0-0. 6-0.0 
n 1 q+0. 06+0. 19+0. 00 
u - iM -0. 01-0. 13-0. 00 


n n«+ - 04 +o-o 1 

u - uo -0.03-0.01 

n nnq+° o 82 +o ooo 
u.uuo_ . 000-0. 000 

n n i +0.04+0.00 
<J.U^_o.o2-0.00 

n nm +0.095+0.000 
u.uui_ .000-0.000 

n cQ+2.02+0.12 

U.O»_o.59_ 01 o 

7 q+7.7+1.6 
' -°-5.4-1.3 


fi Qft+0.03+0.00 
U.UU_ .02-0.00 

n 1 1+0.24+0.01 
u - iy: -0.11-0.01 

n nnv+o on+o- 001 
U.ULW _ .002-0.001 

n 9Q+0.34+0.02 
u - zu -0. 17-0.02 

n Q 4 +l- 8 9+016 
u -"-0.97-0.13 

+5.2+0.7 
*-°-3.6-0.6 


B s ° - K*°r] 
B° I^V 


n 9C+0. 15+0. 49+0. 15+0. 57 
u - zo -0. 13-0. 22-0. 05-0. 15 
n 90+0. 04+0. 46+0. 23+0. 29 
u - zo - 0.04-0. 24-0. 10-0. 15 


q 1 7+0.04+0.10+0.03 
U-l 1 -0.04-0.06-0.01 

n nQ+0-0 2 +0-0 3 +0- 01 
u - U9 -o. 02-0. 02-0. 01 


1 7+0.3+0.2 
L - ' -0.3-0.1 
fi4 +0.33+0.11 
u - ua -0.26-0.11 


n fiO+O- 14+0.07 
u - uz -0. 14-0.08 

n 07+0.35+0.10 

u -°'-0.32-0.08 



Considering the width differences of the two mass eigenstates Bh and Bl, the decay amplitudes squared at time t of 
the state that was a pure B° state at time t = can be parameterized by: 



,-rt 



\A f (t)\* = \(f\B(t))\ 2 = —- (\A f \' + \A f \*) 



cosh 



/ATt 



Hf sinh 



/Am 



V 2 



"J 



+C/ cos(Amt) — Sf sin(Ami) 



(70) 



where Am = ran — ttil > and Ar = Th — is the difference of decay widths for the heavier and lighter B mass 
eigenstates. The time-dependent decay amplitudes squared of another channel B° — > / is obtained from the above 
expression by flipping the signs of the cos(Amt) and sin(Ami) terms. For decays to the CP-conjugate final state, one 
replaces / by /. 

Time-dependent decay amplitudes squared can be simplified in two kinds of cases. In B°-B° system, the small 
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TABLE VI: Direct CP asymmetries (in %) in the B s — > PV decays: the first solution (This work 1) and the second solution 
(This work 2). In both solutions, the chiraly enhanced penguin has been taken into account in B — > VP decay amplitudes. 
The first kinds of uncertainties are from uncertainties in charming penguins and gluonic form factors which are discussed in the 
text; the second kinds of uncertainties are from those in the CKM matrix elements. We also cite theoretical results evaluated 
in QCDF [if]] and PQCD [5^ to make a comparison. 



Modes 


QCDF 


PQCD 


This work 1 


This work 2 


B° K+K- 
B a — > K + K 
B° s -> K°K*° 
B° s -» K*°K° 
B° e - ttV 


o +0.6+8.4+5.1+68.6 
z - z -0. 7-8. 0-5. 9-71.0 
r, h + 1.0+3.8+1.6+47.5 
l - 1.1 -2. 6- 1.3-45.0 
-, 7 +0.4+0.6+0.5+1.4 
-0.5-0.5-0.4-0.8 
n 9+O.O+O.2+O.I+O.2 
u - z -0. 1-0. 3-0. 1-0.1 
97 9+6.1+9.8+2.7+32.0 
z ' -6. 8-5. 6-2. 4-37.1 


or. ,,+2.3+2.8+1.3 
°°- D -2.3-3.5-1.2 

r r +4 4+8 5+5 1 
O0.0_ 4 g _ 9 g _ 2 5 




I Q q + 0.3 + 2.1 + 1.5 

lo - J -0.4-1.7-0.7 


1 1 9+19.1 + 1.3 
A1 ""-16.2— 1.3 
7 i +11.2+0.7 
' • -12.4-0.7 






1 9 q + 11.4+0.8 
iz " 5 -ll. 3-0.8 
q ,,+ 13.0+0.7 
"•"-13.5-0.9 






B° s - p% 
B° - pV 

B s — > U)T) 

B s -» cur] 

B s — ► 

B° s -> <M 


97 0+6.4+9.1+2.6+25.9 
zi -°-6. 7-5. 7-2. 2-28. 4 

90 q+6. 1+10. 3+1. 5+24. 8 
zo - M -7.5- 6.3-1.8-27.5 

4+2.0+30.1+14.6+36.3 
— 2. 1 — 71.2— 44.7— 59.7 

R9 9+15.9+132.3+80.8+122.4 
DZ.^_ 10 2 _ 84.2-46.8- 49.9 


Q 9+1.0+2.8+0.4 
a - z -0.4-2.7-0.7 
9 c o+l. 3+2. 8+3.4 
*SO.O_2, 0-3.6-1. 5 

-1 fi 7+5.8+15.4+0.8 
iD - ' -3.2-19.1-1.7 
7 7+0.4+4.5+9.4 
' • ' -0.1-4.2-0.4 

_-| o+0. 0+0. 6+0.1 
— 0.1—0.6 — 0.2 

7 o+l. 5+1. 2+0.1 
' -"-0. 5-8. 6-0. 4 




U 


Ol q+53.5+2.5 
— 83.2 — 2.6 

a .,+5.3+0.6 
*'*-7.1— 0.6 






o 

u 


Ifi q+W-S+i.e 

— 18.3—1.6 

7 g+5.0+0.8 
' •°-4.9— 0.8 


B t l^n-K*+ 
B° s ->n°K* 

B° 8 ^p°K 
B° s -> K°ui 
B° s - K°</> 


n R +0.2+1.4+0.1+19.9 
u -°-0. 1-1. 7-0. 1-20.1 

ir 7+14.3+13.0+28.4+80.0 
' -16.0-11.6-28.0-59.7 
1 r+0. 4+1. 2+0. 2+12.1 
J -°-0.4-1.4-0.3-12.1 

0/1 7+7.1+14.0+22.8+51.3 
' -5.2-12.4-17.7-52.3 

,iq n+13. 6+18. 0+30. 6+57.7 
^■J- 9 - 13. 4- 18. 2-30. 2-49. 3 
ln o+3. 0+4. 7+3. 7+5.0 
lU.O_2.4_3, 0-4.1-7.5 


1Q O+2.5+2.7+0.9 
1J - u -2.6-3.4-1.4 

A7 -1+7.4+35.5+2.9 
^'•-'-8.7-29.8-7.0 
■>» 9+2.4+2.3+1.2 
J-^-^-2.2-1. 6-0.7 

70 yi+6.4+16.2+2.2 
,y_47 g_3 ,9 

r i +3.2+22.7+3.2 
°*'- L -0.0-15.1-2.0 




Q q+17.2+0.9 
J - J -16.7-0.7 
99 q+33. 1+2.1 
zz - y -40.2-1.9 

11 q + 17.5+1.2 
AA -°-20. 0-1.1 
19 o+20. 1 + 1.0 

94 .+33.7+2.2 
^•1-41. 4-2.0 

q O+5.3+0.3 
J - u -4.7-0.3 


19 , + 17.5+1.1 
iz -^-15.3-1.2 

1 q 4+18.6+0.8 
J.O.fl_ lg _g_ 1 2 

1 n O+9.4+0.9 
1U.O_ 10 2-1.0 

q 9 r+30.7+2.7 
O^.J_23.4_2.9 

1 o 9+16.4+1.2 
AO - Z -17.0-1.7 
9 9+3.0+0.1 
z - z -2. 9-0.1 


B° - K*° V 


/Ifl 9+17.0+24.6+ 7.8+65.9 
11. 5-30. 8- 14. 0-96. 3 
ro c+16. 9+41. 4+19. 9+44.9 
00 -°-ll. 9-11. 7-13. 9-35. 7 


ri 9+6.2+14.1+2.0 
dl - z -6.4-12.4-3.3 
c-1 1+4.6+15.0+3.2 
di - i -6. 6-18. 2-4.1 


9_ „+23. 4+2.0 
zd - ' -22.0-1.3 
qc 9+63.3+3.1 
od - z -49.4-3.8 


fi9 7+28.1+2.6 
OZ. / _22.5-3.9 

q 9 i +22.8+2.6 
oz ' A -23.2-1.7 



width difference Ar can be safely neglected. Thus the first two terms cosh an d smn f^p) * n can ^ c 

reduced to 1 and and the decay amplitudes squared becomes: 



\A f (t)\ 2 = \(f\B(t))f 



(\A f \ 2 + \A f \ 2 ) 1 + C f cos(Amt) - S f sin(Ami) 



In the following, we use the phase convention CP\B°) = \B°) and define the following amplitudes ratios: 



\ - qAf 



\ _<i A f 



and and p are the mixing parameters between B° and B°. The definitions for Cf and Sf are given by: 



C f = 



l-\\ f \ 2 \Atf-\Aff 



1 + IA/l 2 \Af\* + \Aj\* 



Sf 



\A f \* - \Af\z 



1 + IA 



5 / = 2 



im (A/) 
1 + IA/l 2 ' 
Im(A/) 



1 + IA 



/I 



(71) 



(72) 



(73) 



The system of four decay modes defines five asymmetry parameters, Cf, Sf, Cf, Sjr together with the global charge 
asymmetry related to the overall normalization: 



A C p 



\A f \* + \Af\ 2 -\A f \ 2 -\Af\ 2 
\Af\* + \Atf + \Af\* + \Af-r 



(74) 
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TABLE VII: Mixing-induced CP asymmetries in B — » tt* p T decay processes: the first solution (This work 1) and the second 
solution (This work 2). In both cases, the chiraly enhanced penguin has been taken into account. The first kinds of uncertainties 
are from uncertainties in charming penguins which are discussed in the text; the second kinds of uncertainties are from those 
in the CKM matrix elements. We also cite theoretical results evaluated in QCDF approach [h| to make a comparison. 



Parameter 


Exp. 


QCDF 


This work 1 


This work 2 


A C p 

C 

S 
AC 
AS 


-0.13 ±0.04 
0.01 ±0.07 
0.01 ± 0.09 
0.37 ±0.08 

-0.04 ±0.10 


n nn+ - 00 + - 01 + - 00 + - 10 

u.uu_ 00 _o .01-0 .00-0 .10 

n nn+o. oo+o. oi+o. oo+o. 02 

u - uu -0. 00-0. 01-0. 00-0. 02 
n i o+0. 60+0. 04+0. 02+0. 02 
U-lO_0. 65-0. 03-0. 01-0. 01 

n i K +0. 06+0. 23+0. 01+0. 01 
u - lu -0. 07-0. 26-0. 02-0. 02 

n n9 +0. 01+0. 00+0. 00+0. 01 
u - uz -0. 00-0. 01 -0.00-0. 01 


n 19+0.04+0.04 
u ' lz -0.05-0.03 

n n1 +0.13+0.00 
u ' u± -0.12-0.00 
n i i +0.07+0.08 
-0.08-0.13 

n i i +0.12+0.01 
u - ±± -0.13-0.01 
n A 7 +0.08+0.05 
u -^' -0.06-0.04 


n 91 +0.03+0.02 
u - zl -0.02-0.03 

n n1 +0.09+0.00 
u ' u± -0.10-0.00 
n n1 +0.06+0.08 
u - u± -0.07-0.14 

n 19+0.09+0.01 
u - ±z -0.10-0.01 
n ^n+0.05+0.03 
u - q: 'J-0.07-0.03 



One can also use the parameters C = \(C f ± Cf), S = §(S/ + Sf), AC = \(Cf - Cj), AS = \(S f - Sj). If there 
is no direct CP violation, only two independent decay amplitudes squared are left. Thus Acp = 0, Cj = —Cj and 
Sf = —Si which also implies C — and S = 0. If we recall that the CP invariance conditions at the decay amplitudes 
level are Af = Aj and Aj = Af, one can study the following two parameters: 



l%l 2 + l^/l 2 ' ff l^/l 2 + l^/l 2 ' 1 j 



Sometimes, they are considered as more physically intuitive parameters since they characterize direct CP violations. 
In B° — ► p ± 7r T decays, (choosing / = p + n~ and / = p~n + ), we use A+~ parameterizes the direct CP violation 
in decays in which the produced p meson does not contain the spectator quark, while A~^~ parameterizes the direct 
CP violation in decays in which it does. Of course, these two parameters are not independent of the other sets of 
parameters given above, and can be written by: 

= a cp + C ff + A CP AC ff = A C p + C fT + AgpACff 

p7r " l + AC ff - + A CP C ff - ' »* ' -l + AC f f + AcpC fr ' 

Predictions on these parameters are given in table IVIII Most of them are consistent with the data except AC and 
AS*. 

If the final state / is a CP eigenstate, there are only two different amplitudes since |/) = ±|/) and the time- 
dependent decay amplitudes squared can also be simplified. Restricting the final state / to have definite CP-parity, 
the time-dependent decay width for the B — > / decay is: 



r(B u (i) -> /) = e- rt T(B — > /) 



/Art\ /Art 

cosh ^— j+Hf smh ^ — 

cos(Amt) - S f sm(Amt)] . (77) 



The time dependent decay width T(B(t) — > f) is obtained from the above expression by flipping the signs of the 
cos(Ato^) and sin(Ami) terms. In the Bd system, the width differences arc small which can be safely neglected, but 
in the B s system, we expect a much larger decay width difference (AT /T)b s - This is estimated within the standard 



model to have a value (Ar/r) Bs = -0.147 ± 0.060 (6l|, while experimentally (Ar/r) Bs = -0.33+£5? |39|, so that 
both Sf and Hf, can be extracted from the time dependent decays of B s mesons. The definition of the various 
quantities in the above equation are as follows: 

2Im[A] 2Rc[A] 

^-TTpAp f ~ TT|A|2' (78) 
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TABLE VIII: Mixing-induced CP asymmetries Sf in B — > VP decay processes:the first solution (This work 1) and the 
second solution (This work 2). In both cases, the chiraly enhanced penguin has been taken into account. The first kinds 
of uncertainties are from uncertainties in charming penguins and gluonic form factors which are discussed in the text; the 
second kinds of uncertainties are from those in the CKM matrix elements. We also quote the experimental results to make a 
comparison. 



Channel 


Exp. 


This work 1 


This work 2 


B° - P°K S 
B° -> toKs 
B° -» <t>K s 


0-61±o;24 ± °- 09 ± °- 08 
0.48 ±0.24 

0.39 ±0.17 


n or+o.o4+o.oi 

U.OO_ 05 _o 01 

n K1 +0.05+0.02 
u - J1 -0.06-0.02 

0.69 


n cc+o. 02+0.01 

u - du -0.03-0.01 

n sn+°- 02 + 01 

U.OU_ 02-0.OI 

0.69 


B° -> K*-iv + -> Kstt~tv + 

B ' — > K 7T -» K SIT TV 

B° -> A"*°77 -> Astt ^ 
B° -> A*V ^ Ksir°T)' 




n Q n+0.04+0.01 
U.90_q .07-0 .02 

/-\ pt-*-^0 04^0 02 
u - JZ -0.05-0.02 
n 7 p+0.01+0.01 
u - ' J -0.01-0.01 

n 7C+0.07+0.01 
u - '"-0.06—0.01 


n O.+0.06+0.03 
u - oy: -0.07-0.03 

u - ' a -0.02-0.01 
n R ,+0.01+0.00 
u -°*-0.01-0.00 

n fifi+o.o4+o.oo 

U.OO_ .05-0.00 


S(B° 7T°p°) 

S(B° -> 7r°o;) 


0.12 ±0.38 


n i 1+0.14+0.10 
u - iJ --0.14-0.15 
fl O7+0.44+0.02 
u -°' -0.00-0.01 


n iQ+0.14+0.10 
u - iy -0.14-0.15 
n 7O+0.36+0.07 
u - ' z -1.54-0.11 


B° -> A 
B° -> pV 
S° -> a;? 7 




fl oc+0. 15+0.03 
U.OO_2.03-0.07 

n 7Q +0.20+0.05 
u - '^-1.73-0.09 
n i O+0.19+0.10 
u - iz -0.20-0.17 
n no+0.59+0.10 
u - z,3 -1.10-0.10 


n 9Q +0.36+0.09 
u - zy -0.44-0.15 

QQ+0.22+0.09 
U.OO_ 1 .24-0.14 

n R +0.14+0.10 
u - iD -0.15-0.15 
n 9 7+0.17+0.09 
u - z ' -0.33-0.14 



with 

q A(B^f) 

x -^- p A(B^Ty (79) 

where rjf is +1(— 1) for a CP-even (CP-odd) final state /. q/p = e~ 2l P for the Bd system while q/p = e +2le for the B s 
system where e = arg[— V c bV ts V* s V t * b ]. With the convention arg[y c h] = argfvy = 0, the parameter can be reduced to 
e = arg[— VtsV^]. For b — > s transition induced B° decays, the ratios of decay amplitudes j^^/j are almost real and 
thus Sf ~ sin(2/3). These channels provide a good way to measure sin(2/3). Experimentalists often use the following 
parameters in b — > s transitions: 

-'V 5 '-" 2 l + |AP ' -IM— 2 1 + IAP ' (80) 
while the latter parameter is only defined for the i3° — B® system. Although the K*° meson is not a CP eigenstate, 
its daughter-mesons Ksir behave as CP eigenstates. Thus we also give the predictions on mixing-induced CP 
asymmetries in the decays involving a K*° meson and other related decays. Results for these parameters are collected 
in table IVIlTl and tabic HXl where predictions on decays with branching ratios smaller than 10~ 7 are omitted. 

After studying the two simplified cases, we come to the time-dependent CP asymmetries in B® — > K* + K~, where 
the final state are not CP eigenstate and the width difference of B® — B® can not be neglected cither. In the following, 
we choose / = K* + K~ and / = K + K*~ . One needs to consider two additional CP asymmetries: 

_ Re(A / ) _ 2 ifc(Aj) (m 

which can be redefined as: H = Hi ^ H f and AH = f 2 1 . Our predictions for these parameters are given in table IXl 
but we have not considered the global charge asymmetries because of the presence of Ar. These predictions will be 
tested at the forthcoming LHCb experiments 
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TABLE IX: Mixing-induced CP asymmetries (S/)s 3 and (#/)b s in B s — * PV decays. Results obtained in the PQCD 
approach [ss| are also collected here; the errors for these entries correspond to the uncertainties in the input hadronic quantities 
(charming penguins and the two form factors £ s and C,j a ), and the CKM matrix elements, respectively. 



Modes 



PQCD 



This work 1 



This work 2 



B° - A 



B a -> tvrj 



+0.01+0.08+0.02 
0.01-0.09-0.03 

3+0.00+0.01+0.01 
'-0.00-0.03-0.00 
+0.06+0.14+0.01 
-0.06-0.16-0.01 



-0.07 



0.15 



n Qo+o. oi+o.oi+o. oo 

u - yo -0. 01-0. 03-0. 00 



n oq+o.oo+0.04 

u - oa -0.00-0.05 

_n A CL+0.00+0.09 
u -* J -0.00-0.10 



n Qn+ ' 0(J + ' 02 

U.»U_ 00 _o .03 

n 4^+0.00+0.05 
u - y:y: -0.00-0.05 



1.00 



+0.00+0.00 
-0.06-0.01 



0.60 



+0.30+0.03 
-0.53-0.03 



_q -j^g+0. 00+0. 10+0. 04 

0.95 



0.00-0.12-0.05 
+0.01+0.01+0.01 
-0.00-0.02-0.02 



-n nzL+°' 41 + ' 08 

u - uy: -0.32-0.08 
q gg+0.00+0.02 

0.32 



1.60-0.02 
+0.67+0.06 
-1.29-0.06 



0.02 



+0.01+0.02+0.00 
-0.03-0.08-0.00 



0.99 



+0.01+0.01+0.00 



-0.11 
0.99 



-0.01-0.06-0.00 
+0.01+0.04+0.02 
0.00-0.04-0.03 
+0.00+0.00+0.00 
0.00-0.00-0.00 



n fi9 +0.41+0.08 
u - u ^-0.18-0.12 



+0.16+0.11 
-0.20-0.06 



-0.79 
-0.25 



n on+°' 20 +0-02 

u -° u -0.36-0.02 

-0 4i+°- 75 + ' 10 
-0.91 

q g3+0.C 

-0.37 



0.75-0.15 
+0.82+0.08 
-0.08-0.04 



-1.23+0.10 



-0.03 
1.00 
0.00 
1.00 



+0.02+0.07+0.01 
-0.01-0.20-0.02 

+0.00+0.00+0.00 
-0.00-0.01-0.00 
+0.00+0.02+0.00 
-0.00-0.02-0.00 
+0.00+0.00+0.00 
-0.00-0.00-0.02 



n q-7+2. 12+0.05 
u - y '-0.00-0.02 
+0.43+0.04 
-0.15-0.04 



-0.39 



-0.98-0.04 
7+1.37+0.09 
-0.65-0.10 
1-0.04+0.01 
-0.00-0.00 
n nn+0-32+0.12 
U.U»_ .22-0.08 

+0.35+0.02 



-1.00 



n Qn+ 014 + 02 

U.»U_ .24-0.02 

-0 07 +0,06+0,01 

l.ooi 



0.06-0.01 

0.00+0.00 
0.01-0.00 



0.23 
0.96 
0.10 
0.99 



+0.04+0.01 



+0.07+0.01 



+0.01+0.00 
-0.01-0.00 



b: -» i<< 



B° s ^p°Ks 



B s -> K s uj 



-0.72 
-0.69 
-0 57+ - 22 +°- 51 + - 02 

-0.36 
-0.63 
-0.57 



0.17-0.39-0.05 
+0.10+0.46+0.04 
0.13-0.15-0.04 

+0.09+0.28+0.01 
0.09-0.11-0.02 
+0.11+0.31+0.02 
-0.13-0.38-0.02 



q gg+0.04+0.01 
-1.00 

99+ C 
0.04 



0.03-0.01 
+0.00+0.00 
-0.00-0.00 



-0 i3+°- 02 + - 01 
-0.99 



0.05-0.01 
+0.13+0.09 
-0.11-0.13 
_ n i 1+0.28+0.18 
u - 11 -0.22-0.14 
n QC+0-02+0.01 
u -»°-0. 16-0.03 



-0.03 



0.02-0.01 

+0.00+0.00 
-0.00-0.00 

+0.22+0.17 
-0.17-0.12 
+0.05+0.01 
-0.13-0.02 
+0.02+0.00 



0.95 
0.98 



n n7 +0.11+0.08 
u - u '-0.09-0.12 



B" -> K*+n- -»• K s n + TT- 



DO . IV^*0^0 , jy _0 u 
ti a — > K 7T — > KsTT 7T 



B° -► K*°r) -> K s n°r) 



B° -» K*"r/ -► K s n° V ' 



0.98 



+0.01+0.01 
-0.04-0.02 



n iR+o.ii+0.09 

u - lu -0.09-0.13 

-0 07 +0 ' 26+0 18 
0.97 



0.22-0.14 
+0.01+0.01 
-0.15-0.02 
+0.06+0.02 
-0.09-0.03 
n oo +0. 15+0.08 
u - zz -0.14-0.11 



0.94 



-0.94 



+0.33+0.03 



0.01 



+0.45+0.16 



+0.11+0.15 
u -°°-0.09-0.11 

QQ+0 03+0.03 
u - ao -0.07-0.07 
Q4+0.03+0.02 
u - a ^-0.05-0.04 

_ n nn+ ' 12 + ' 07 
U - JU -0.09-0.10 

_fl 77+0- 2 3+0.04 
u - ' ' -0.16-0.03 

n 1n +0. 26+0.11 
U.1U_ .22-0.11 

n 7O+0.15+0.04 
u - ' z -0. 16-0.05 
_ n R n+0.20+0.05 
u - uz -0. 16-0.06 



G. Isospin asymmetries and U-spin asymmetries 

Currently, there are many experimental methods to measure CKM angles: a, (3 and 7. But in order to reduce the 
uncertainties, a good way is to use SU(3) symmetry although this will induce the errors from SU(3) symmetry breaking 
effect. Here we will present some tests on this kind of symmetry breaking, although the flavor SU(3) symmetry for 
B — > P ', B — > V form factors and various charming penguins are used. 
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TABLE X: Mixing-induced CP asymmetries in — > K* + K~ decay processes: the first solution (This work 1) and the second 
solution (This work 2). In both predictions, we have included the chiraly enhanced penguin and chosen / = K* + K~ . The 
first kinds of uncertainties are from uncertainties in charming penguins which are discussed in the text; the second kinds of 
uncertainties are from those in the CKM matrix elements. 



Parameter 


This work 1 


This work 2 


C 

S 
H 
AC 
AS 
AH 


n n9 +o. 10+0.00 
u - uz -0. 11-0.00 
n n9 +0.07+0.01 
u - uz -0.07-0.01 
n Q9 +0.02+0.02 
u - MZ -0.04-0.02 

-fl nq +0 - 11+0 - 01 

u - ua -0. 10-0.01 
n QQ+0.07+0.04 
u -°°-0.07-0.04 

n m +0.04+0.00 
u - ui -0.02-0.00 


n m +0.09+0.00 
u - ul -0.09-0.00 
n no +0.05+0.01 
u - uz -0.05-0.00 
n Q1 +0.02+0.02 
U.»l_o.03-0.02 
n i +0.09+0.01 
U.ll_o.09-0.01 
n .1+0.05+0.03 
u - 4:1 -0.05-0.03 

n nl +0.02+0.00 
u - ui -0.02-0.00 



TABLE XI: Two kinds of results for the ratios Ris in B — > 7T7t and B — > %p decays, together with the predictions in QCDF [lOj] 
and experimental data evaluated using the results of branching fractions. The first kinds of uncertainties are from uncertainties 
in charming penguins as discussed in the text; the second kinds of uncertainties are from those in the CKM matrix elements. 





Exp. 


QCDF 


This work 1 


This work 2 


Rx 
R 2 
R 3 
R4 
R5 


2 RQ+ - 54 
z - oy -0.53 

9 91 +0.37 
z ' zl -0.37 

J--OO_0.46 
OR+0-80 

u -°'-0.33 


9 q Q +0. 31+0. 04+0. 15+0.05 
Z.Oy_ .25-0.08-0.12-0.11 

9 nfi +0. 40+0. 53+0. 12+0.03 
z ' uo -0. 30-0. 36-0. 09-0. 06 
-1 no+0. 18+0. 82+0. 03+0. 02 
J-JO-0. 17-0.59-0.04-0.05 
n Ao+0. 04+0. 15+0. 45+0. 23 
u - 4:Z -0. 04-0. 11 -0.21 -0.20 
n 99+0.07+0.08+0.23+0.14 
u - zz -0. 08-0. OC-0. 12-0. 12 


-■ 09+0.15+0.10 

■■■•"■"-0.12— 0.12 
i -17+O.I3+O.O6 
L - 1 ' -0.11-0.07 

1 90+0. 12+0. 08 
J-^S-O.IO-O.IO 
9 oQ+0.21+0.02 
_.OO_ 20-0 .02 

1 91+O.O5+O.O2 
i - ZJ --0.05-0.03 


1 fi/l +0.22+0.05 
J-°^-0.19-0.06 
1 t-9+0.17+0.07 
■■-■"^-0. 15-0.09 
1 e^+0.07+0.11 
J- dq; -0.09-0.07 
1 9O+0.04+0.02 
i ' ZJ -0.04-0.02 
1 nQ +0.08+0.04 
J- UJ -0.08-0.03 



In the B — > 7T7t and B — > irp system, one often uses the following ratios 10( : 



Ri 


_ r(B° - 


-* 7T + p~ 


) 


R 2 = 


r{B° - 


■+ W + p~) +T(B° -> 7r-p+) 


~ T(B° - 


-> 7T+7T 






2r(S° -> 7T+7T-) 


R 3 


_ T(B° - 
~ T{B° - 


+ 7T + p" 


) 


i?4 = 


2T(B- 
T(B° - 


^7T-p°) _ 2r(B-^7rV) 1 
->n-p+) ' ° r(£?°^7r+p-) 



(82) 

where the partial decay widths are CP averaged. Our predictions are given in table IXI| where we have used the 
experimental results on branching ratios to evaluate the ratios and these values are collected as experimental results. 
The predictions in QCDF approach arc also collected in this table. In B° — > tt + it~ and B° — ^ ir + p~ . tree operators 
dominate. If we only consider the tree operators, R\ becomes ratios of decay constants: R\ = (f p / f^) 2 ~ 2. Our 
predictions are smaller than 2 for both solutions. In the first solution, the ratio is much smaller which is mainly caused 
by charming penguin terms: A p f gives a constructive contribution to the decay width of B° — > 7r + 7r~ while AY C P gives 
a destructive contribution to F(B° — > 7r + p~). In the second solution, the deviation of R\ from 2 is not too large as the 
phase of A pp is almost the same as A v c p . Ri and R b are larger than predictions in QCDF approach and the present 
experimental data. B~ — ► 7r~p° contains two different contributions from tree operators: color-allowed contribution 
with p~ emitted; color-suppressed contribution with 7r _ emitted. In QCDF approach, the second contribution is small 
and the first contribution is related to tree operators in B~ — > n~p + . Neglecting the color-suppressed contribution 
and contributions from penguin operators, R4 is equal to zero. In SCET, color-suppressed tree operators can give 
sizable contributions as we have discussed. Thus the branching ratio of B~ — ► ir~ p° is enhanced which can give a 
large value for i? 4 . The analysis is also similar for the ratio R5. 

In — > K*~ir + , B® — > K + p~ , B% — > K~p + and B® — > K* + ir~ , the branching ratios are very different from each 
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other due to the differing strong and weak phases entering in the tree and penguin amplitudes. However, as shown 
by Gronau 62j , the two relevant products of the CKM matrix elements entering in the expressions for the direct CP 
asymmetries in these decays are equal, and, as stressed by Lipkin subsequently, the final states in these decays 
are charge conjugates, and the strong interactions being charge-conjugation invariant, the direct CP asymmetry in 
B® — > K + -k~ can be related to the well-measured CP asymmetry in the decay B® — ► A~7r + using U-spin symmetry. 



In this symmetry limit, we have 62|, |63[: 



\A(B° S -> t:+K*-)\ 2 ~ \A(B° S - n-K*+)\ 2 = \A(B d -> p+R-)\ 2 - \A(B d - p~K + )\ 2 , (83) 

CPV P ; CPK s ' BR{B a d -> r(S s ) v ; 

Following the suggestions in the literature, we can test these equations and search for possible new physics effects 
which would likely violate these relations. Accordingly, one can define the following parameters: 

\A(B S -> n+K*-)\ 2 - \A(B S -> ^A*+)| 2 = BK(B S -> 7r-A*+)Ag^ -> 7r-A'*+)r(B d ) 
6 ~ |^l(i? rf ^p-A+)| 2 -|A( J B d ^p+A-)| 2 BK{B^K-p+)A% r P (B^K-p+)T(B s ) ' 1 j 

-> p+A-) i?A(g s ^7r+A*-) Tjgrf) 
1 n-K*+) + BR(B d p+K-) r(B s ) ' 1 j 

= \A(B S -> p+A-)| 2 - |A(g a -> p-A+)| 2 _ ff7£(ff s -> p-A'+)Ag^(B a -> p~ A+)r(£? d ) 
7 " ' \A(B d -+ p-K+)\ 2 - \A{B d p+A-)| 2 6^(5 A*-tt+)^(B - K*-p+)r(B s ) ' 1 j 



Ag^(B d ->7r+A*-) BR(B a ->p-K+) r(B d ) 
A(4 p (B s -> p-K+) BR(B d -f 7T+A*-) ' r(S s ) ' 



We also consider B° — > 7r + p , B . — > A+A* , B° — > 7r p + and B° — > K* + K which are related by U-spin 
transformation and define the following ratios: 

. |A(fl a -> A*+A-)| 2 - |A(B a -> A*-A+)| 2 _ g7e(g s -> A+A*-)Ag£(g a -> K+K*-)r(B d ) 

8 " l^^rf^p+TT-)! 2 -!^^^^)! 2 -> 7T+p-)Ag- (B - 7T+p-)r(5 s ) ' [ ' 



A d j r P (B d -> p-7T+) BR(B S ^K*-K+) r(B d ) 



A$ P (B S -> A*-A+) BR(B d -> p-ir+) t(B s ) ' 



(90) 



|A(B a -> A+A*~)| 2 - |A(fl a -> A~A*+)| 2 _ g^(B s -> A~A*+)Ag^(B a -> A~A*+)T(B d ) 
9 " |^(B d ^+p-)| 2 -|A(£ d ^ 7 r- /9 +)| 2 B7i(B-.7r-p+)i4*J(S-.7r-p+)r(B a ) ' 1 ' 

A% P (B d ^n-p+) BR(B S ^K~K*+) r{B d ) 
4 A d j p (B s A"A*+) + BR(B d ^n-p+) t{B s ) ' 1 j 

In the flavor SU(3) symmetry limit, the ratios are i? = — 1 and A is zero. Using the first solution for the 16 inputs, 
we obtain the following values: 

_ n ns+o-02+o.oi 

u - uo -0.04-0.01' 

_ n m +0.00+0.00 

u - ul -0.01-0.00' 

n i i +0.06+0.03 
u - li -0.05-0.02; 

o OO+0.12+0.06 /qo\ 
u -°°-0.11-0.05' V yo / 



R 6 = -0.89, 


Ax = 


R 7 = -0.99, 


A 2 = 


A 8 = -1.11, 


A 3 = 


R 9 = -1.24, 


A 4 = 


are omitted here. 


Our 


R 6 = -0.87, 


Ai = 


i? 7 = -0.99, 


A 2 = 


A 8 = -1.10, 


A 3 = 


i? 9 = -1.25, 


A 4 = 



n in +0.03+0.02 
-U.±U_ 05-0 .02> 

n ni +0.00+0.00 
-U.UJ._ .00-0.00' 

nQ +0.03+0.01 
■ ua -0.02-0.01' 



0.33i°;i?l ;°ol. (94) 
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Since the form factors and charming penguins are assumed to the respect flavor SU(3) symmetry, the small deviations 
for the ratios R and A are reasonable. 



V. COMPARISONS WITH THE PQCD APPROACH 




(a) 



d{s) 



4 



4 



(b) 



FIG. 2: Feynman diagrams for the (S — P)(S + P) annihilation operators in PQCD approach and charming penguins in SCET. 



PQCD approach is based on kx factorization, where one keeps the intrinsic transverse momentum of quark degrees 
of freedom. The intrinsic transverse momentum can smear the end-point singularities which often appear in collinear 
factorization. Rcsummation of double logarithms results in the Sudakov factor which suppresses contributions from 
the end-point region to make the PQCD approach more self-consistent. This approach can explain many problems to 
achieve great successes. Currently, radiative corrections 5lU64l.l65l.l66l| and power corrections in 1/mj, 67], [68[ in this 



approach are under studies. In PQCD approach, annihilation diagrams can be directly calculated. Among them, the 
(S — P)(S + P) annihilation penguin operators (from the Ficrz transformation of (V — A)(V + A) operators) are the 
most important one. According to the power counting in PQCD approach, annihilation diagrams arc suppressed by 
J^-QCD/fnb but the suppression for (S—P)(S+P) annihilation penguin operators is 2r x . This factor is comparable with 
1. Thus annihilations play a very important role in PQCD approach. Phenomenologically, the large annihilations can 



tt+tt- and B° -> K--k + 



69| , the polarization 



explain the correct branching ratios and direct CP asymmetries of B° 

problem of B — > 4>K* |to| . etc. In Fig.[2ja), we draw the Feynman diagrams for this term. Comparing with charming 
penguins, we can see they have the same topologies in flavor space. So generally speaking, charming penguins in 
SCET as shown in Fig.[2fb) have the same role with (S — P)(S + P) annihilation penguin operators in PQCD. Both 
of them are essential to explain the branching ratios in these two different approaches. But there are indeed some 
differences in predictions on other parameters such as direct CP asymmetries and mixing-induced CP asymmetries. 

First of all, the CKM matrix elements associated with charming penguins and (S — P)(S + P) annihilation penguin 
operators are different. If we consider B decays in which a b quark annihilates, the (S — P)(S+P) annihilation penguin 
operators are proportional to VtbV t * Dl while charming penguins are proportional to V c bV* D . The differences in the 
CKM matrix elements will affect direct CP asymmetries and mixing-induced CP asymmetries sizably. For example, in 
B° — > (f)Ks decay, the mixing-induced CP asymmetries in SCET are dramatically different from predictions in PQCD 
approach. In the SCET framework, there is no contributions from tree operators to B s — > (f>Ks at tree level and 
penguin operators are much smaller than charming penguins. As the CKM matrix element V c bV* D for the charming 
penguin is real, the parameter A defined in Eq. (|79|) becomes A = — e +2le , where we have neglected contributions from 
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penguin operators. Thus in SCET the two parameters Sf and Hf are given by: 

S f = - sin(2e) = -0.03, H f = - cos(2e) = -1.00. (95) 

In PQCD approach, the CKM matrix clement for the (S — P)(S + P) annihilation penguin operators is VtbV^ d which 
gives A = - e +2^+2i/3. 

Sf = - sin(2e + 2/3) = -0.72, H f = - cos(2e + 2(3) = -0.69. (96) 

The differences in the mixing-induced CP asymmetries between SCET and PQCD will be tested at the future exper- 
iments. 

In PQCD approach, contributions from the (S — P)(S + P) annihilation penguin operators can be calculated using 
perturbation theory. These contributions are expressed as the convolution of light-cone distribution amplitudes and 
a hard kernel. We can also include SU(3) symmetry breaking effects in the calculation in PQCD approach. In SCET, 
charming penguins are from the charm quark loops. Since the charm quark is heavy, one can not factorize charming 
penguins (see Ref. , \l\ for another point of view) . Thus charming penguins are non-perturbative in nature 

which is similar with the final state interactions 
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73] . In the present work based on SCET, we have assumed SU(3) 
symmetries for the contributions from charming penguins. The magnitudes and strong phases of charming penguins 
can not be calculated using perturbation theory which obtained by fitting the experimental data. 

The third difference is the magnitudes of charming penguins in SCET and contributions from the (S — P)(S + P) 
annihilation penguin operators in PQCD approach. This difference arises from the different power counting in the 
two approaches. We take b — > s transitions to illustrate the difference. In PQCD approach, the (S — P)(S + P) 
annihilation penguins are enhanced to be of the same order with penguins in emission diagrams. In SCET, charming 
penguins are more important. Comparing the values given in Eq. (|50[) . Eq. (|54[) and Eq. (|60[) . we can see charming 
penguins in SCET always larger than contributions from emission penguin diagrams. 

In PQCD approach, the (S — P)(S + P) annihilation penguin operators are chiraly enhanced and the dominant 
contribution is from the imaginary part. The main strong phases in PQCD approach which are essential to explain the 
large CP asymmetries in many channels are also produced through from these operators. But in SCET, as we have 
shown in Eq. (T5D|) and Eq. ([5"^]) . strong phases of charming penguins are not too large. Accordingly, our predictions 
on direct CP asymmetries are small compared with predictions in PQCD approach. 



VI. CONCLUSIONS 



We provide the analysis of charmless two-body B — > VP decays under the framework of soft-collincar-cffcctivc 
theory. Besides the leading power contributions, we also take some power corrections (chiraly enhanced penguins) 
into account. In the present framework, decay amplitudes of B — > PP and B — > VP decay channels can be expressed 
as functions of 16 non-perturbative inputs: 6 form factors and 5 complex (10 real) charming penguins. Using the 
B — > PP and B — > VP experimental data on branching fractions and CP asymmetry variables, we find two kinds 
of solutions in % 2 fit for these 16 non-perturbative inputs. Chiraly enhanced penguin could change some charming 
penguins sizably, since they have the same topology with each other. However, most of other non-perturbative inputs 
and predictions on branching ratios and CP asymmetries are not changed too much. With the two sets of inputs, we 
predict branching fractions and CP asymmetries. Agreements and differences with results in QCD factorization and 
perturbativc QCD approach arc also analyzed. Our conclusions are as follows: 
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• In color- allowed processes such as B° — > n^p^ decays, tree operators provide the dominant contributions. Our 
predictions on branching fractions are smaller than the ones calculated in the QCDF approach and PQCD 
approach. The main reason is that: both B — > P and B — > V form factors in SCET arc smaller. B° — > ir°p° 
and other color-suppressed channels are predicted with a larger branching ratios in SCET, because the hard 
scattering form factors (j' V are comparable with £ p,v who also have a large Wilson coefficients. The large 
branching ratios for B — > ir°p° are consistent with the experimental data. 

• b — > s decay processes such as B — > wK*, B — ► pif and the corresponding £? s decays are dominated by 
contributions from charming penguins. Since we have assumed flavor SU(3) symmetry for charming penguins, 
branching fractions of b — > s transition decays can be estimated by analyzing the corresponding charming 
penguin terms. Decays with iso-singlet mesons rj and rj' are slightly different since there exists cancelations 
between different charming penguins. 

• In the PQCD approach, annihilation diagrams do not suffer from the endpoint singularity problem, which can 
be directly calculated. Among the three kinds of penguin operators, the (S — P)(S + P) operators are most 
important which provide the main strong phase in the PQCD approach. In the SCET framework, charming 
penguins play an important role especially in b — ► s transitions. The (S — P)(S + P) annihilations have the 
same topology with charming penguin. Besides the commons, there exists many differences in these two objects 
including weak phases, magnitudes, strong phases, SU(3) symmetry property and factorization property. These 
differences will mainly affect the direct CP asymmetries and time-dependent CP asymmetry variables. 
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APPENDIX A: EXPRESSIONS FOR HARD KERNELS 

For explicit decay channels, the hard kernels depend the Lorentz structure and flavor structures. They can be 
evaluated using the Wilson coefficients given in Eq. (fT6|) and Eq. p7|) . In this appendix, we intend to write the decay 
amplitudes in a compact form. In doing it, the following meson matrices are required: 



B- = (1,0,0), B° = (0,1,0), B° s = (0,0,1), 









c 


0^ 












[oo 


0^ 










(o 





(A 






M n+ = 


M p+ = 


1 


c 





• 


M K+ 


= M K ,+ = 












M K o 


= M K »o 













! 










c 


°J 
(l 





o\ 








[l 


Q J 


(l 


0^ 








1 


°) 

(o 





o\ 


V2M^o 


= V2M p o = 







-1 







V2M Vq = 









1 




= M+ 





















1° 





°J 












V° 


J 










\ {) 







M n - = 


M p - = Mj + 




M K - 


= M K , 






M R0 


= M R , 


= M£o, 















32 



we also need the following matrices: 





(l 









M 




M 


S u = 











, A d = 


i 


, A s = 







v° 





V 








w 



(A2) 



Using the meson matrices, one can write the hard kernels appearing in B — > M1M2 decays as: 



Ti = c f 1 BM 2 S u M 1 A f + {4 ±4)BAf 2 A / Tr[5 u A/i] 

+c{BM 2 M 1 A f + {c{ ± c^BMaA-^TrfMr], 
Ti 9 = cfs^A/iA^rlMal + ^ic^BA/Tr^MxlTr^] 

+c{BAf 1 A / Tr[il/ 2 ] + (c{ ± 4)BA. f Tt[M 1 ]Tc[M 2 ], 
Tf = clBM 2 K f Ti[M x ], T? B = c? g BlJTSc\Mi]'Tr[Mzl 

Tu = TM b{), T 1Jg = T X M b{), Tfj = Tf(c{ -> &{), Tf Jg = Tf s (c{ -> ftf). 



(A3) 



If the emitted meson M 2 is a pseudoscalar, — C3 and — Cg in are used. But for vector meson emission, we use 
plus signs in the combinations. 

Using meson matrices, the charming penguins responsible for B — > M±M 2 decays can be determined in the same 
way. If the charming penguins in B — > PP decays are considered, the master equation is: 



a m iM2 = BM2MlA I ' A pp + BM 1 A fr Yi[M 2 }A 



pp 

ccg J 



(A4) 



where the A ccg term is only responsible for the iso-singlet mesons rj q and rj s . In B — > VP decays, the charming 
penguins are: 



I A/i AI 2 



BM 2 M 1 A f A V J + BMiM 2 A f A P J + BMiA f Tr[M 2 ]A 



VP 

ccg ! 



(A5) 



where we take Afi as a vector meson and M 2 as a pseudo-scalar meson. 

The master equations for hard kernels for chiraly enhanced penguins are given by: 
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